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ABSTRACT 

We present calculations of magnetic potential functions associated with the perturbation of 
Saturn's planetary magnetic field by a rotating, equatorially-situated disc of plasma. Such 
structures are central to the dynamics of the rapidly rotating magnetospheres of Saturn and 
Jupiter. They are 'fed' internally by sources of plasma from moons such as Enceladus (Sat- 
urn) and lo (Jupiter). For thesemodels, we use a scaled form of Caudal's Euler potentials for 
the Jovian magnetodisc field (Caudal 1986). In this formalism, the magnetic field is assumed 
to be azimuthally symmetric about the planet's axis of rotation, and plasma temperature is 
constant along a field line. We perturb the dipole potential ( 'homogeneous' solution) by using 
simplified distributions of plasma pressure and angular velocity for both planets, based on 
observations by the Cassini (Saturn) and Voyager (Jupiter) spacecraft. Our results quantify 
the degree of radial 'stretching' exerted on the dipolar field lines through the plasma's rota- 
tional motion and pressure. A simplified version of the field model, the 'homogeneous disc', 
can be used to easily estimate the distance of transition in the outer magnetosphere between 
pressure-dominated and centrifugally-dominateddisc structure. We comment on the degree of 
equatorial confinement as represented by the scale height associated with disc ions of varying 
mass and temperature. For the case of Saturn, we identify the principal forces which contribute 
to the magnetodisc current and make comparisons between the field structure predicted by the 
model and magnetic field measurements from the Cassini spacecraft. For the case of Jupiter, 
we reproduce Caudal's original calculation in order to validate our model implementation. 
We also show that compared to Saturn, where plasma pressure gradient is, on average, weaker 
than centrifugal force, the outer plasmadisc of Jupiter is clearly a pressure-dominated struc- 
ture. 

Key words: (magnetohydrodynamics) MHD — plasmas — methods: numerical — planets 
and satellites: general 



1 INTRODUCTION 

Jupiter and Saturn are not only the largest pl anets in our Sola r sys- 
tem, they are also the most rapid rotators. iGledhilll ( Il967h first 
pointed out the important consequences of these properties for 
Jupiter's magnetosphere. The rotational period of the planet is ap- 
proximately 10 h, and as a result the gravitational (Fg) and cen- 
trifugal (Fc) forces associated with corotating plasma in Jupiter's 
magnetosphere are equal at an equatorial distance of ~ 2.3 Rj 
from the planet's centre (here, we denote Jupiter's radius as Rj ~ 
71000 km). At the orbit of lo, situated at 6Rj, centrifugal force 
exceeds gravitational by a factor of nearly 20. Saturn's radius 
(Rs ~ 60000 km) and rotational period (~ 10.75 h) modify these 
distances to 1.7 Rg {Fg = Fc) and 4.7 Rs {Fg ^ Fc/20), the latter 
being 0.7 Rs outside the orbital radius of the icy moon Enceladus. 
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Evidently, centrifugal force is an important factor for determining 
the structure of the outer magnetospheres of these planets. 

IGledhilll jl967h showed that the action of centrifugal force in 
Jupiter's rapidly-rotating magnetosphere tends to confine magne- 
tospheric plasma towards the equatorial plane, where the planet's 
assumed dipolar field lines reach their maximum radial distances. 
A disc-like magnetospheric structure was thus anticipated, and in- 
deed observed by the first spacecraft to visit the Jovian system. 
Pioneers 10 and 11 ( Smith etal. 1974, 1975). The near-equatorial 
magnetic field structure seen by these spacecraft was very differ- 
ent from that of a rotating dipole throughout the region referred to 
as the middle magnetosphere, situated at distances ~ 20-50 R.j. 
The magnetometer observations in this region showed a periodic 
pattern of largely radial field direction alternating with intervals 
having a north-south (meridionally-directed) field. These data were 
interpreted as periodic encounters with a rotating, disc-like cur- 
rent sheet. The highly radial field in this picture is a signature of 
magnetic field lines resembling a dipole pattern that has been radi- 
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ally 'stretched' outwards near the magnetic equatorial plane. Such 
a magnetic geometry would be associated with an inward Lorentz 
force ~ J0 X B, part of which is required to provide the cen- 
tripetal acceleration for the rotating magnetospheric plasma (here 
denotes azimuthal current density and B magnetic field). 

Observations by Galileo ( lKivelsonet"aDll997h confirmed the 
persistence of this magnetodisc structure, and examined its re- 
sponse to changing solar wind conditions. During the late inbound 
pass of the Galileo insertion orbit, the magnetic field measurements 
indicated th at a strong compres sion of the magnetosphere had 
taken place (*Kivels on et alj[l997l) . This compression witnessed by 
Galileo resulted in an increase of the meridional field Be by a factor 
of ~ 2 when compared with the data from the Pioneer 10 outbound 
segment over the middle magnetospheric region at 30-50 Rj . Like 
Galileo, Pioneer 10 outbound was a near-equatorial swathe situ- 
ated at a local time near dawn. Unlike Galileo, however, the Pi- 
oneer 10 observations, acquired 22 years earlier, indicated a rela- 
tively quiescent magnetosphere. The conclusion was that the mag- 
netospheric compression at the time of the Galileo insertion had 
caused an increase in meridional field Be by squeezing the mag- 
netic flux threading the magnetospheric plasma into a smaller vol- 
ume (the change in location of the dawn magnetopause was in- 
ferred to be 40 Rj inward). In addition, the periodic field signa- 
tures in meridional and radial field (Be and Br) seen by Galileo 
indicated a thicker plasma sheet within the magnetodisc structure, 
as one would expect for a strongly-compressed magnetosphere. 

The an alogous behav i our for the magnetodisc at Saturn was 
explored by I Arridge et alj ( l2008bl) , who took advantage of many 
orbits of magnetometer measurements from the Cassini spacecraft 
in order to investigate the relationship between magnetosphere size 
(as represented by the subsolar magnetopause standoff distance 
Rmp) and the degree to which the radial field Br dominated the 
field measurements seen in the outer magnetosphere. For the ob- 
servations considered in this study, Cassini was typically situated 
on near-equatorial orbits just outside, and south of, the current 
sheet. This analysis revealed that, under conditions of low solar 
wind dynamic pressure (corresponding to Rmp > 23 Rs), the 
magnetic field due to Saturn's ring current (i.e. azimuthal current) 
dominates the planetary internal field in the outer magnetosphere, 
with the combination of the two producing the magnetodisc struc- 
ture. For a more compressed magnetosphere, however, the dayside 
field became strongly dipolar, with magnetodisc geometry surviv- 
ing only on the magnetosphere 's nightside and flanks. The Kro- 
nian magnetodisc may thus essentially disappear on the dayside 
under appropriate conditions, and is therefore even more sensitive 
than Ju piter's magnetodisc in r esponse to upstream solar wind con- 
ditions jArridge et al]|2008bl) . ISunce et alj \200 §) arrived at simi- 
lar conclusions by modelling the response of the ring current and 
its magnetic moment for different magnetospheric configurations. 
This was done by the application of an empirical field model (CAN 
disc) to different orbits of Cassini magnetometer data. We describe 
this field model in further detail below. 

Connernev et al. ( 1981, 1981) provided the first detailed mod- 
elling of the ring current which supports the magnetodisc field at 
Saturn, and applied this model to the magnetometer data from the 
Voyager spacecraft encounters. The magnetic field in the Conner- 
ney, or CAN, model is computed by assuming, a priori, an az- 
imuthally symmetric distribution of current which is confined to 
an annular disc of uniform thickness, extending between inner and 
outer edges at specified radial distances. Current density in this 
model is assumed to be inversely proportional to cylindrical ra- 
dial distance (J^ oc 1/p). This formalism has been employed in 



several studies of the structure of Saturn's ring current, usually 
based on fitting in situ magnetic field measurements from space - 
craft (e.g. Conne rnev et"^ l ll983h : lGiampieri & Dougherty! ( l2004h : 
[Bunce et al. ( 200^ )). 

The study bv lSunce et alj j2007h emphasised that the current 
which flows in the magnetodisc current sheet is a macroscopic man- 
ifestation of the microscopic drift motions of charged particles in 
the plasma. These authors examined the contribution of two types 
of azimuthal particle drift to the magnetic moment of the ring cur- 
rent: (i) the magnetic gradient drift exhibited by particles of finite 
thermal energy whose guiding centre moves in response to changes 
in field strength experienced during individual gyrations, (ii) the 
inertial drift associated with the centrifugal force in a frame which 
corotates with the local plasma flow. They showed that, for typi- 
cal magnetospheric conditions at Saturn, the heavier (water-group) 
ions may generate a much stronger inertial current at distances be- 
yond ~ 10 Rs due to their rotational kinetic energy exceeding typ- 
ical thermal energy. 

Theoretical and empirical magnetic field models for the ring 
current at Jupiter a nd Sa t urn proposed by var i ous authors (e.g. 
'Gleeson & Axford| ( Il976h : lOoertz et al.l j 19761) : IConnemev etaU 
(1981, 1983)) have proved to be valuable tools for determining the 
global length scales and inte nsity of t he cur rent which supports the 
magnetodisc field structure. ICaudall i fT98^ pointed out that the a 
priori current distributions used in such models cannot be used 
to infer, unambiguously, the dynamical properties of the plasma 
in which the current flows. In particular, determining the relative 
importance of the plasma pressure gradient and centrifugal forces 
in generating the plasma current and magnetodisc field, requires a 
different approach which incorporates a knowledge of the plasma 
prop erties. 

ICauda]| ( ll986h developed a formalism in which Jupiter's mag- 
netic field structure was modelled by solving a magnetostatic equa- 
tion representing dynamical equilibrium, i.e. a uniformly zero vec- 
tor sum for all of the aforementioned forces throughout a speci- 
fied region. This solution was then used to infer the global dis- 
tribution of current which was consistent with the equatorial dis- 
tribution of plasma properties such as angu lar velocity, tempera- 
ture, density and composition. Caudal (1986) used observations by 
the Voyager spacecraft (Con nernev et al. 1981; Bagenal & Sulliyag 
1 198 ll : iKrimigis et al.ll 1981^ to constrain this equatorial plasma in- 
formation, which, in his formalism, acts as a boundary condition for 
inferring the global plasma properties. The resulting current distri- 
bution from such a calculation has a more realistic global structure 
than the uniformly thick, annular disc used in the empirical models. 
By including the effects of both plasma therm al pressure and cen- 
trifugal force in his formalism, ICaudall [m^ naturally extended 
previous investigations of the distortion of the planetary magnetic 
field which assumed a cold plasma with negligible th ermal energ y 
compared to rotational kinetic energy (e.g. Hin&Carb'^jl978l) ). 

The main purpose of this paper is to adapt the formalism by 
ICaudall jl98^ in order to model the magnetodisc of Saturn. For 
the required equatorial plasma properties, we use t he latest ob- 
serva tions by Cassini Plasma Spectrometer (CAPS) jYoung et al.l 
|2004|) and Cassini M agnetospheric Imaging Instrument (MIMI) 
I Krimigis et al.ll20(j4h . The framework, assumptions and inputs for 
the model are summarised in fj2] For the sake of completeness, we 
provi de a derivation of the magnetostatic solution cited by Cauda| 
( Il986i) in Appendix |A] This derivation is not published elsewhere, 
to the best of our knowledge. Its inclusion here serves as a starting 
point for discussion of a toy model for the plasmadisc described 
in ^2.2\ This model has a very simplified structure in terms of its 
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plasma properties, but serves as a useful illustration of the compe- 
tition between plasma pressure and centrifugal forces in determin- 
ing magnetodisc structure. Detailed magnetodisc models for Sat- 
urn are presented in !j3]and compared with Cassini magnetometer 
(MAG) data from e quatorial and high-latitud e orbits. A description 
of MAG is given in lDoughertv et al] ( l2004h . These model outputs 
are also compared with those of the best-fitting CAN discs. We 
conclude with a summary and discussion in ij?) 



2 MODEL FRAMEWORK 

2.1 Magnetic Field Geometry and Force Balance 

We adopt the formalism of lCaudaj ( Il986h . and express the magnetic 
field components associated with an axially symmetric plasma dis- 
tribution as gradients of a magnetic Euler potential a. The value of 
a is constant along any magnetic field line. It is also constant over 
any axisymmetric shell of field lines (flux shell). The change in a 
between flux shells is simply related to the magnetic flux contained 
between them (i.e. it is a flux function, see ij2.3.2| l. With this as- 
sumption, the magnetic field radial component Br and meridional 
component Bg are 



Br = 



Ba = 



da 



1 da 
r sin 6 dr ' 



(1) 



where 9 denotes colatitude with respect to the planetary rotation 
axis (assumed coincident with the magnetic axis), and r is radial 
distance from planet centre (in units of planetary radii). The unit 
of a in our 'normalised' system is Boa, the product of the equa- 
torial magnetic field Bo at the planet surface and the planet radius 
a. The adopted values and corresponding scales for relevant phys- 
ical quantities at both Jupiter and Saturn are shown in Table[T]and 
Table IB 1 1 (Appendix iBt . Unless otherwise stated, we shall use this 
dimensionless form of Caudal's original equations in order to eas- 
ily compare the degree to which different plasm a discs may distor t 
the internal field of their parent planets (see also lVasvliunasI ( 1200 Sl) 

and 

Caudall ( Il986l) examined the condition of general force bal- 
ance in the rotating plasma 



J X B = VP — nrriiUj pep, 



(2) 



where J is current density, B is magnetic field and p — r sin 8 is 
cylindrical radial distance from the axis (ep is the corresponding 
unit vector). Plasma properties are pressure P, temperature T (as- 
sumed isotropic and constant along field lines), ion number density 
n, mean ion mass rrii and angular velocity cj. This equation repre- 
sents balance between magnetic force on the left side, and pressure 
gradient plus centrifugal force on the right. We have not included 
the minor contribution to plasma mass from the electrons, but do 
include their contribution to plasma pressure. 

Caudal ( 1986) used the definitions of field and current ex- 
pressed as functions of a. When these expressions are substituted 
into the force balance condition (Eq. the result is the following 
partial differential equation 

d^a d^a 

where p is the cosine of colatitude i.e. p = cos 0. 

The 'source function' g is determined by the global distribu- 
tion of plasma pressure and angular velocity. ICaudall ( ll986l) pointed 



out that g could be used to derive the azimuthal current density 
according to 



J4>{r,p) = 



r sin p 



(4) 



Force balance in the direction parallel to the magnetic field im- 
plies that the global values of these quantities are derivable from 
their equatorial values and the shape of the magnetic field lines. 
This is why there is a general dependence of g upon a. Ca udall 
(1986) derived an analytical expression which could be used to cal- 
culate the solution for a. We have included a full derivation of this 
expression in Appendix |A] The form involves the use of Jacobi 
polynomials, which occur as solutions of the homogeneous version 
of Eq. An important solution in this class is the dipole poten- 
tial Qdip ~ (1 ^ A'^)/''- In practice, we start with a pure dipole 
potential and then 'perturb' it using Caudal's iterative method: at 
every iteration, the solution a„ is used to evaluate g and thus the 
'next' solution an+i . We stopped our calculations when the differ- 
ence between successive iterations was at most 0.5 per cent. We 
describe the various inputs used for our Saturn model calculations 
in i]2.3l These are based on a variety of observational studies em- 
ploying data taken by the Cassini spacecraft. Before investigating 
these Caudalian disc models for Saturn, we shall examine a sim- 
ple toy model which may be used to predict the effect of a rotating 
plasma disc upon the zeroth-order (largest-scale) perturbation to 
the dipole potential. 



2.2 Toy Model for a Planetary Magnetodisc 

We begin our investigation by examining a very simplified model of 
disc structure. In this model, we assume that the disc has a uniform 
plasma P parameter denoted Ph, associated with the thermal energy 
of a hot population. We also assume the presence of an isothermal 
cold population containing most of the plasma mass, but a negli- 
gibly small fraction of the total pressure, with uniform plasma /3 
denoted /3c. Proc eeding under t hese assumptions, it is straightfor- 
ward to show that lCaudallfl98^ 's expression for the plasma source 
function may be written as 



g{r,p,a) = p^ 



da 



+ p exp 



2P 



PcO 



(5) 



where Ph and Pc denote hot and cold plasma pressure, and the sub- 
script is used to refer to the quantity evaluated at the equatorial 
crossing point of the magnetic field line, i.e. the magnetically con- 
jugate point for which p — 0. The relation between source func- 
tion and current density Eq. l|4) allows the identification of the two 
terms on the right-hand side of Eq. l|5j as quantities proportional 
to the individual contributions to total current density which arise 
from hot pres sure gradient, and from centrifugal force. Following 
ICaudall 119861) , we assume that the hot component exhibits uniform 
pressure Ph all the way along a given field line, while the cold 
component's pressure is concentrated towards the equatorial plane, 
according to the exponential factor in Eq. l|5j- The symbol £ thus 
represents a scale length a ssociated with the cylindrical radial dis- 
tance p, and is defined by ( ICaudaliri986h 

f - (6) 

for a quasi-neutral plasma containing singly charged ions and elec- 
trons, a represents planetary radius and is used here to transform to 
our normalised system (see Table [TJ. 

It is worth emphasising here that this expression for £ arises 
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Table 1. Physical units used in the normalised (dimensionless) system for both planets. 



Planet 


Radius (a) 


Magnetic Field (Bo) 


Magnetic Flux (Boa^) 


Pressure (Bq/^q) 


Angular Velocity (cjq) 


Saturn 


60280 km 


21160 nT 


77 GWb 


0.00036 Pa 


27r/10.78 radh^l 


Jupiter 


71492 km 


428000 nT 


2187 GWb 


0.146 Pa 


27r/9.925 radh^l 



from lCaudall ( Il986h 's mathematical treatment of general force bal- 
ance, but it is also a natural consequence of field-aligned force bal- 
ance for both ions and electrons in a quasi-neutral plasma. Ions 
and electrons in this formalism are implicitly subject to a simpli- 
fied ambipolar electric field which acts to distribute both types of 
particle with an equal scale length £ given by Eq. l|6j. We can see 
how this arises by considering the equations for field-aligned force 
balance, for both ions and electrons, which take into account par- 
ticle pressure, centrifugal force (projected along the magnetic field 
direction) and the presence of the ambipolar electric potential $||. 
This derivation of the scale length £ is summarised in Appendix Icl 
Using this definition of plasma scale length, it can be shown 
that the exponential factor appearing in Eq. l|5j has an argument 
which contains the ratio of an ion's kinetic energy of rotation to its 
thermal energy. Thus our hot plasma component is defined by ion 
thermal energies which are large compared to the kinetic energy of 
rotation at angular velocity lu and a consequent scale length which 
is effectively infinite (large compared to magnetospheric flux tube 
length). On the other hand, the cold plasma component contains 
ions with much smaller thermal energies, which cannot compete as 
effectively with the centrifugal potential in maintaining plasma at 
locations high above the equatorial plane. We shall see in the later 
sections that the typical scale lengths for the cold plasma at Saturn 
are a few planetary radii, a distance which is small compared to the 
flux tube lengths in the planet's outer magnetosphere. 

If we now make the assumption for the global magnetic field 
Beo = (where e.g. x = 3 for a dipolar geometry), then by def- 
inition the dependence of the normalised magnetic pressure along 
the equator is given by ^Pq'^^- It follows that P^q = ^PhPo'^^ 
and Pco ~ ^PcPo^^. Replacing the operator with the equiva- 
and assuming uniform T and lo, we may transform 



lent 

the expression for the plasma source function from Eq. Q into the 



following form: 

gir,p.,o:) 



2P 



exp 



pI 



2^2 



(7) 



It is straightforward to show, using Eq. l|6]l, that the term 
{j3cp'^ /2l'^) is equivalent to the ratio of rotational kinetic energy 
density to magnetic pressure. It may therefore be thought of as a 
plasma 'beta' for bulk rotation, rather than random ion motions. If 
we consider the equatorial location of any given flux tube (p = po), 
then we find that the hot and cold plasma contributions to the source 
function are equal at an equatorial radial distance pr given by: 



pl^2x£^{fiH/Pa). 



(8) 



Beyond the transition distance pr the rotational kinetic en- 
ergy of the plasma exceeds its thermal energy. Therefore, centrifu- 
gal force (in a corotating frame) dominates pressure gradients for 
distances p ^ pr in determining both the magnetospheric cur- 
rent; and the distortion to the planetary dipole field required to 
maintain the magnetodisc's dynamic equilibrium. Conversely, for 
p <S pT, the role of rotation is less important and plasma pres- 
sure determines disc structure. Interestingly, pr may conceivably 
exceed the standoff distance of the dayside magnetopause under 



conditions where: (i) hot plasma (3 is very high compared to the 
cold plasma, (ii) plasma angular velocity is adequately low, or (iii) 
for a given temperature of cold plasma, its density is small (such 
that the quantity £'^ / fic becomes very large). We shall see in the 
following sections that average magnetospheric conditions at Sat- 
urn may yield transition distances at or inside the magnetopause 
(outer magnetosphere), while for Jupiter the magnetospheric cur- 
rent arises predominantly from hot pressure gradients. 

We shall now investigate the relatively simple expression for 
the magnetic potential of the homogeneous disc model. We make 
use of the following equality, valid for dipolar magnetic field lines: 

po=r/(l-p'). (9) 

Eq. ([9} indicates that dipolar field lines have parabolic shapes in 
the (r, p) co-ordinate plane. By definition, p^ = r^(l — /i^) and 
thus the argument of the exponential factor in ^ may be expressed 
in terms of r and p. The full representation of the plasma source 
function in these co-ordinates can be derived as (eliminating po 
using Eq. 



gD(r,p) = r ^(1 - p ) 



2NX+1 



2^2 



exp 



2^2 (1_^2)2 



(10) 

The source function for our simplified disc structure, Eq. jlOh 
has been derived using a dipolar magnetic field. Therefore it is an 
appropriate form for magnetic fields which may be decomposed 
into the planetary dipole plus a small perturbation (in the sense that 
the perturbation is everywhere small compared to the dipole field 
strength). In order to provide an approximation for this perturbation 
field due to the plasma disc, we shall calculate only the zeroth-order 
terms in the expansion for the magnetodisc potential described by 
Eq. l lA26t . By 'zeroth-order', we mean the terms involving the or- 
thonormal basis function (1 - p^) Pq '^(p) = (V3/2)(l - p^), 
which involves the Jacobi polynomial of order zero. This function 
represents the largest angular scale of the magnetodisc potential 
(p^ dependence). The planetary dipole field has this angular de- 
pendence and is thus included in the zeroth-order solution. 

If we use Eq. l lA26t to calculate the zeroth-order part of the 
potential, in conjunction with the explicit form of the source func- 
tion for our homogeneous disc. Eq. ilO\ . the resulting expression 



ao(r,p) = (1 - p )/r 



1+ u go 



3 f°° -1 \ 
{u)du + r / u go{u)du\ 



(11) 

where we have defined go, the zeroth-order coefficient of the source 
function go, in terms of radial distance: 



1 

9o{r) = - J gD{r,p)dp. 



(12) 



The first factor in Eq. i ll U is the unperturbed dipole poten- 
tial, and the integral terms in the second factor (enclosed by square 
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Figure 1. Geometry of magnetic field lines for the zeroth-order homoge- 
neous disc models (see text). The solid black lines indicate contours which 
are spaced by uniform intervals in the logarithm of the magnetic potential. 
The dark grey and light grey lines indicate regions which lie between the 
same values of magnetic potential in each panel. 

brackets) represent the lowest-order (largest angular scale) pertur- 
bations due to the presence of the model plasma disc. It is evident 
from the integral limits that any location which lies outside the disc 
plasma will still experience a magnetic field due to all of the remote 
disc currents flowing within the radial distance of such a point. 

Fig. [T] shows contours of equal magnetic potential (i.e. mag- 
netic field lines) for the function qq evaluated for three examples 
of the homogeneous disc model. The first model is a hot disc with 
Ph = ^ and no cold population {J3c = 0); the second is a cold disc 
with Ph = G and 13c — 0.2; and the third is a combined disc with 
~ 0.5 and /3c — 0.1. All the homogeneous discs are assumed to 
be in perfect corotation with the planet, have uniform length scale 
£ — 1 and uniform field strength index x ~ 3. The disc models 
have had their structures truncated by setting their source function 
to zero for regions which are magnetically conjugate to equatorial 
distances po < 5 and po > 35. The field lines of a vacuum dipole 
are also shown for comparison in the top panel. The closely-spaced 
dark grey and light grey field lines indicate regions in each panel 
which cover the same intervals in magnetic potential. This may also 
be considered as an interval in magnetic flux, since the total mag- 
netic flux <I> B threading the region from cylindrical radial distance 
p to infinity is a simple multiple of magnetic potential evaluated at 



this distance: ^b{p) = 2-Ka{p, Z) (where the vertical co-ordinate 
Z — r cos = rp). This relation follo ws from the de finition of the 
magnetic field components in Eq. lUll ( ICaudaliri986l) . We see that 
the presence of the disc tends to 'inflate' the dipolar magnetic field 
lines and shift them to larger equatorial crossing distances. For the 
flux tubes highlighted, we see that this effect is quite pronounced 
for the dipolar field lines which cross 10-12 planet radii (Rp) at 
the equator (shaded light grey). These flux tubes are displaced to 
equatorial distances 16-25 Rp by the disc models. By contrast, the 
flux tubes shaded dark grey, situated near ~ 5 Rp, do not undergo 
as great a distortion in the presence of the disc. This part of the 
field may be considered as a 'rigid', inner magnetosphere domi- 
nated by the internal planetary field. We also note that the cold disc 
model produces outer magnetospheric field lines which are notice- 
ably more oblate in shape compared to the other discs: a field line 
crossing at a given equatorial distance does not rise as far above the 
equator in the cold disc model. This property reflects the tendency 
of the cold plasma disc to be concentrated near the equatorial plane 
(according to the scale length £), thus producing stronger distortion 
in the near-equatorial segments of the model field. 

The outer magnetospheric flux tubes in the present example 
(field lines shaded light grey in Fig.[T) also become spread out over 
a larger radial distance compared to the dipole configuration. In 
the upper panel of Fig. |2] we compare the corresponding equato- 
rial profile of field strength between the disc models and the dipole 
field. For all disc models, the ratio of total to dipole field strength 
increases monotonically with distance, indicating that the magne- 
todisc field is more uniform than the dipole. We also see that each 
model has a characteristic distance which separates an inner region 
where the field strength ratio B/Bdip < 1 from an outer region 
where this ratio exceeds unity. This feature is observed in actual 
planetary ring currents, and arises from the finite extent of the cur- 
rent region and the solenoid-like nature of the corresponding mag- 
netic field (e.g. Sozou & Windld (md ) and Figures 1, 4 and 5 of 
Arridge et al. (2008b)). Near the inner boundary of the equatorial 
disc current, the vertical magnetic field generated by this current 
alone opposes the planetary field, while the opposite is true near the 
outer disc boundary, where the disc field enhances that of the plan- 
etary dipole. We thus expect and find that the ratio of total magnetic 
field to that of a pure dipole monotonically increases from values 
less than unity near the inner edge of the disc (~ 5Rp) to values 
larger than unity near the outer edge (~ 35 Rp). 

We now consider the radial profiles of the main forces in- 
volved in the dynamic equilibrium of the homogeneous disc struc- 
ture. We do not expect these forces to be in perfect balance for 
our zeroth-order model, since this is only one component of the 
many required to retrieve the full solution, and corresponds to the 
largest angular scales of the problem (sin^ 6 dependence). We plot 
the equatorial centrifugal force, magnetic forces and plasma pres- 
sure gradients as a function of distance in the lower panel of Fig.|2] 
for the homogeneous, combined disc model described above. The 
total magnetic force (J x B) is the sum of the magnetic pressure 
gradient and the curvature force. When these two components have 
equal magnitude and opposite sign, the total magnetic force is zero. 
In the figure, the absolute value of the negative curvature force is 
displayed. Thus it is the vertical difference between this curve and 
that for the magnetic pressure gradient which indicates the mag- 
nitude of the total magnetic force. The plots show us that, for the 
outer region where the ratio B/Bdip > 1 and p > 12 Rp, the 
force of highest magnitude is that due to magnetic curvature, fol- 
lowed by centrifugal force (factor of < 2 smaller than curvature 
force) and magnetic pressure gradient (factor of > 5 smaller than 
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Figure 2. Upper panel: field strength in the equatorial plane, relative to 
that of a pure dipole field, as a function of radial distance p, for various 
homogeneous disc models (see legend and text). Lower panel: forces per 
unit volume (dimensionless) in the equatorial plane of the homogeneous, 
combined disc model, labelled according to line style. 



curvature force). The hot plasma pressure has a gradient about half 
the magnitude or less of that for magnetic pressure, while the weak- 
est force in this model is the gradient due to cold plasma pressure. 
We see that the sum total of the forces is less than one percent 
of the curvature force at p = 12 Rp, and less than ten percent 
for p — 10-16 Rp. This aspect of the total force profile is an in- 
dication of the degree to which dynamic equilibrium is maintained 
within the zeroth-order component of the full solution. 

The transition distance for the combined disc model is pT ~ 
5.5 Rp, using Eq. ([8}. It is also evident from Fig. [2] that, for dis- 
tances much greater than this value, we are by definition in the 
part of the magnetosphere where centrifugal force dominates hot 
plasma pressure, and where the force balance principally involves 
the centrifugal and magnetic curvature forces. As we approach 
the transition distance from the outer disc, the other forces due 
to plasma pressure and magnetic pressure become comparable to 
curvature force; and thus become more significant in determining 
stress balance and disc structure. 



mean mass rrii between these two limits, given by: 

nwinw + npmp 



nw + np 



(13) 



where the symbol n denotes number density, with obvious sub- 
scripts indicating water group and proton components. 

In order to capture the behaviour of rri i as a function of ra - 
dial distance, we employed the formulae of IWilson et al.l ( |2008|) . 
who determined and fitted density moment s for water group io ns 
and protons using observations by CAPS. I Wilson et al.l (|2pOS|)'s 
observations sampled five equatorial spacecraft orbits in the dis- 
tance range ~ 5.5-11 Rs between October 2005 and April 2006. 
The orbits were chosen as mission segments during which CAPS 
ion mass spectrometer obtained sufficient coverage of the plasma 
particle distribution to allow reliable computation of moments. For 
the purposes of our modelling, w e used the Ga ussian fits to water 
ion density and proton density bv lWilson et al.l l|2008) to compute 
the following number fraction of protons as a function of radial 
distance: 



fM{p) 



nw + riH 1 + {Aw/AH)exp[{BH - Bw)p'^] 

V- 15' 



/a/(p)=0.1 



1 — tanh 



+ 0.8, 



(14) 



where the function parameters provided bv lWilson et alj ( 1200 Sl) are 
Aw = 161.5 cm"^, Ah = 8.3 cm"^, Bw = 0.042 Rs"^ and 
Bh = 0.031 Rs~'^- Since these fitted functions are based on ob- 
servations in the distance range 5-12 Rg, we used the hyperbolic 
tangent function /m to place the additional constraint that the pro- 
ton number fraction approaches 80 per cent in the o uter magneto- 
sphere . This plasma composition was determined by I Arridge et al.l 
( l2007t) to provide good agreement with both electron densities ob- 
served by Cassini in the outer magnetosphere and the surface mass 
density of the Kronian plasma disc. The latter quantity was de- 
duced in the same study from the analysis of magnetic signatures 
of transient excursions by the spacecraft into the magnetodisc cur- 
rent sheet. Fig. [3] sh ows plots of the fitted composition profile from 
IWilson et al.l ( |2008|) ; the extrapolation of this profile beyond the 
range of validity (5-12 Rg); the profile used in the current work; 
and the multiplier function /m • 



2.3 Model Inputs and Boundary Conditions 

In this section, we summarise the model inputs we have used to 
determine the equatorial boundary conditions for our calculations 
of the Kronian magnetodisc field. These inputs have been drawn 
from a variety of observational studies using the plasma instru- 
ments aboard the Cassini spacecraft. They fall into the following 
four categories, each described in more detail in the subsections be- 
low: (i) composition, (ii) temperature and density, (iii) hot plasma 
pressure; and (iv) rotation. 



2.3.1 Equatorial Plasma Composition 

The equatorial composition of the cold plasma at Saturn is required 
in the model for the computation of the scale length I (Eq. (|6}), 
assumed constant along magnetic field lines. The composition of 
the plasma is determined by the relative densities of water group 
ions (mas s mw = 18 a mu) and protons (mass mp — 1 amu). 
Following ICaud 3 ( Il986h we represent the disc ions as having a 



2.3.2 Equatorial Plasma Temperature and Density 

For the cold plasma population, the contribution to the source func- 
tion g (Eq. l|3j) takes the form (see also Eq. l|5]l): 

g^ir,p,a)=p e.p^^^ji^— + ^J, (15) 

where the geometry of the magnetic field, represented by a, de- 
termines the mapping between p and po along a field line, as well 
as the equatorial field strength Bgo ■ In order to compute the scale 
length £ in the equatorial plane, we require the equatorial distribu- 
tion of plasma temperature. Strictly speaking the latter should be 
the field-parallel temperature, since (. is associated with force bal- 
ance parallel to the magnetic field. In addition, a knowledge of both 
plasma temperature and density is required to specify the equato- 
rial pressure Pcp. To satisfy the se requirements, we appealed to 
the study by Iwilson et al.l ( |2008|) which provided tabulated mea- 
surements of both parallel and perpendicular temperature for the 
thermal water group ions and protons at Saturn. These tabulations 
contain average temperatures over intervals of radial width 0.5 Rg 
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Cold plasma composition 




Wilson et al. (2008) 
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Figure 3. Profiles of proton number fraction in the cold disc plasma. 
Da shed line: profile de termined from the fits to plasma density moments 
by IWilson et al.l )2008l) . Solid line: profile used in this work. Dotted fine: 
the multiplying function used to constrain the outer magnetospheric com- 
position (se e text). Vertical lines indicate the range of validity of the fitted 
functions of lwilson et all j2008l) . 



near the planet's equatorial plane, extending between radial dis- 
tances 5.5-10 Rs- To obtain total plasma temperatures for mod- 
elling purposes, we began b y combining the tabu lated ion and pro- 
ton temperature values from lWilson et alj ^20081) as follows: 



nwTww + rbHTnw 

nw + nH 
nwTwi^ + tihThi 



nw + riH 
+ 2Tj_ 



(16) 



where nw and nH ar e the respective wat er-group ion and proton 
number densities from lWilson et al. I( l2008h (see d3n i: the symbol 
T represents temperatures; and the subscripts || and _L are associ- 
ated with thermal motions parallel and perpendicular to the mag- 
netic field. The quantities Ty and T± are average parallel and per- 
pendicular temperatures for the cold plasma (weighted by number 
density between protons and water-group ions), while Tc is an ap- 
propriately weighted mean. In Fig. |4] we show plots of the radial 
profiles of fcbTii and ktTc, expre ssed in units of eV. 

While the data provided bv lWilsonet"ai] ( l2008h are valuable 
for our work, we still need to assign temperatures to those regions 
of the magnetosphere outside the reach of this study, i.e. po < 5 Rs 
and po > 10 Rs. In order to do this, we have assumed that the 
individual water group ion and proto n temperature s in th ese re- 
gions are equal to those measured by Iwilson et al.l ( 120081) at the 
closest relevant points (i.e. at 5.5 Rs and 10 Rs respectively). We 
then compute extrapolated total temperatures using the weighted 
sum (according to plasma composition) given by Eq. l ll6t . Even 
though the individual temperatures of the heavy ions and protons 
are assumed fixed in this extrapolation, the variation in the plasma 
composition produces a total temperature which steadily decreases 
with distance beyond 10 Rs. This behaviour is due to the increas- 
ing fraction of the relatively cold protons in the plasma in the more 
distant magnetosphere. We show the extrapolated total parallel and 
mean plasma temperature as dashed lines in Fig. [4] To obtain the 
final, realistic input temperature profiles for the cold plasma, we 
applied second-order polynomial fits to the data-derived profiles of 



and Tc- We believe that suc h an approach is just ified in light of 
the fact that the observations bv lWilsonet"ai] ( l2008h show variabil- 
ity in temperature moments within their 0.5 Rs bins, typically by 
factors between two and five, even for data from the same orbit. The 
final fitted profiles are shown as grey curves in Fig.|4] We note that 
our fitted values for Ty in the range 10-15 Rs are somewhat lower 
than the value of a few hundred e V corresponding to the observa- 
tions of iMcAndrews et ^ ('2009). However, the large-scale trend 
of our fit agrees with these data and our fitted temperatures are of 
similar order of magnitude. We aim to incorporate further plasma 
temperature measurements as they become available. 

We used the profiles of Ty in conjunction with our mean ion 
mass mi ( ^2.3. l\ to compute the plasma scale length. We used the 
Tc profiles in order to compute the cold plasma pressure in the 
equatorial plane, accordi ng to the follo wing formula, adapted to 
dimensionless form from lCaudall ( ll986h : 



Po{a) = 2NL{a){ktTcr/Vw{a). 



(17) 



Here, the dependence upon the local value of a (i.e. the particular 
field line) is indicated for the dimensionless quantities Vw and Nl ; 
these are, respectively, the weighted unit flux tube volume and the 
flux tube content. Considering these two quantities for the moment, 
their definition is based on the usual concept of the unit flux tube 
volume, which we define in our normalised system as: 



V{a) = 



ds/B, 



(18) 



where the integral is taken along a magnetic field line of length 
SB between its southern and northern ionospheric footpoints; ds 
is an element of arc length along the field line; and B is lo- 
cal field strength. Given the relation between equatorial field 
strength and the increment in magnetic potential (see derivation 
of Eq. l[7]l), it follows that 2nV{a) \da\ represents the normalised 
volume between two magnetic shells corresponding to the inter- 
val [a, a + da]. This same volume is threaded by an increment 
27r \da\ of normalised magnetic flux. 

Using this definition of the unit flux tube volume, we can con- 
struct the weighted flux tube volume as follows: 



Vw{a) = 



exp 



2P 



ds/B, 



(19) 



where the exponential weighting factor is a consequence of field- 
aligned pressure balance for the cold rotating plasma (see also 
Eq. (O). The flux tube content Nl is defined as the number of cold 
ions per unit of magnetic flux. That is, the quantity 2nNL{a) \ da\ 
is the number of ions within the volume bounded by the magnetic 
shells corresponding to the interval [a, a + da] in magnetic poten- 
tial. The factor (kbTc)* in Eq. dlVI l is the dimensionless form of the 
thermal energy corresponding to the averaged plasma temperature 
Tc defined in Eq. l ll6t . We obtain this factor through division by 
the energy scaling factor in Table IB 1 1 

By specifying a profile of flux tube content in the Caudalian 
model rather than density, it is more straightforward to mimic re- 
alistic changes associated with a plasma which is 'frozen-in' to 
the magnetospheric field. Our profile for the flux tube content 
wa s obtained by fitting estim ates of this quantity from the work 
by McAn drews et al.l ^20091) . extended to cover inner magneto- 
spheric regions 4-10 Rs (H. J. McAndrews, private communica- 
tion). These authors used Cassini nightside near-equatorial plasma 
observati ons by CAPS in conjun ction with the magnetospheric field 
model bv lKhuranaetalJ j2006l) in order to estimate Nl through 
a force balance relation similar to Eq. dlVl l. We have fitted the 
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Cold Plasma: Equatorial Density Profiles 
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Figure 4. Profiles of pai'allel (top panel) and mean (lower panel) total tem- 
perature for the cold disc plasma. Solid lines with dots: total plasma temper- 
atures derived from the moments f or water group ions and protons, deter- 
mined from Cassini plasma data bv lwilson et alj j20O8h (see text). Dashed 
lines: extrapolated temperature profiles, derived by assuming that the indi- 
vidual ion species have temperatures eq ual to those at the nearest location 
in the tabulation of Iwilson et alj j20O8l) . Solid grey lines: the final model 
inputs, obtained by second-order polynomial fits to the temperature profiles 
derived from the data. 



flux tube content measurements with two Gaussian functions, con- 
strained to meet continuously at 6.5 Rs- The entire fitted profile for 
Nl used in the model is shown in Fig. |5] In order to bring this 
profile of flux tube content into reasonable agre ement with the den- 
sity moments provided by Wilson et al. ( 2008) (who studied orbits 
distinct from those used by Mc Andrews et al. (2009)), we multi- 
plied them by a smooth correction function whose shape, but not 
absolute scale, is also shown in Fig. [5] 



2.3.3 Hot Plasma Pressure 



ICaudall ( [T98^ assumed that the hot magnetospheric plasmas filled 
each flux tube such that each flux tube can be characterised by a 
particular equatorial pressure Pho and volume per unit flux V , re- 
ferred to as the hot plasma approximation. Using the ideal gas equa- 
tion per unit flux one can show that the product Pho V is equal to 
NokbTh where A'^o is the number of ions per unit flux a nd kbTh 
is the mean kinetic energy of the ions (see also ^3^. ICaudall 
used observations from the Jovian magnetosphere to show 
that kbTh did not vary appreciably with L. In the absence of plasma 
sources the time-stationary radial (cross-L) transport of pla sma can 
be des cribed using a one dimensional diffusion equation. ICaudall 
( Il986l and references therein) showed that the plasma tended to 
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Figure 5. Profiles of cold plasma flux tube content (top panel) and num- 
ber density (lower panel) in Saturn's equatorial plane. The top panel shows: 
radial profiles of flux tube content A^^ in ions per Weber f or the Gaussian 
flts to the original observations of lMcAndrews et j2009h (grey line); the 
shape of the multiplying function used to modify this profile (dashed line, 
see text); and the final form of the profile used for the disc models (black 
line). The vertical line shows the division between regions where different 
Gaussian profiles were used to fit the Nj^ data. The lower panel shows the 
agreement between the io n number density profile from the observational 
fits of Wilson et al. with that derived from a disc model with mag- 

netopause radius Rmp = 25 Rs. 



have a uniform distribution in L when the rate of loss of particles 
due to pitch angle scattering into the loss cone was negligible com- 
pared to the rates for cross-L transport. This reasoning led Caud3 
i ll986l) to conclude that PhoV and No{L) were independent of L 
and hence that under rapid radial diffusion the hot plasma in the 
Jovian magnetosphere behaves isother mally rather tha n adiabati- 
cally: PhoV^ = const where 7 = 1. ICaudall (m^ used pub- 
lished energetic particle pressures and magnetic field models in or- 
der to show that the particles did indeed behave isothermally be- 
yond ~ 18 Rj, but adiabati c on smaller L-shells. In further work, 
ICaudal & ConnernevI Il989l) made 7 a free parameter in a fit of the 
model to Voyager magnetometer data. They found that 7 — 0.88 
beyond 9 Rj, suggesting the presence of non-adiabatic cooling pro- 
cesses during inward diffusion, losses, and violations of the first 
and second adiabatic invariants. 

Following Caudal ( 1986), we parametrised the distribution of 
hot plasma pressure in our model by appealing to observations, us- 
ing the same PhoV' = const theore tical framework. Th e data re- 
quired were taken from the study bv lSergis et al.l ( l2007h , who de- 
termined pressure moments for ions with energy > 3keV from 
the measurements of MIMI. The observations presented by these 
authors were acquired within the distance interval 5 < p < 20 Rs 
over eleven consecutive near-equatorial orbits of the spacecraft be- 
tween late 2005 and early 2006. An important result to emerge 
was that, within this 'hot population', particles with energy > 
lOkeV carried half of the total pressure, but contributed only 
~ 10 per cent of the total number density. We shall see in the later 
sections that the hot plasma pressure for typical conditions at Sat- 
urn may exceed that of the colder population (see i]2.3.2t by up 
to an order of magnitude; it was therefore important to include a 
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representation of this liot pressure component in our magnetodisc 
model's source function. 

We used an empirical magnetic field model to determine the 
unit flux tube volume V{a) as a function of p. The empirical model 
comprised an un-tilted dipole and CAN current sheet, where the 
parameters of the model current shee t were dependent o n the dis- 
tance to the subsolar magnetopause teunce et al.]|2007l) . We note 
in passing that the results of this analysis are not significa ntly al - 
tered by using an alternative model, such as iKhurana et all (l2006h. 
The se cond-order fits to plasma /3 as a function of L by Sergis et al] 
( I2OO7 !) were then used to provide values of the hot plasma pressure 
Pho at the equator, using the model magnetic field strength to cal- 
culate pressure from plasma /3. The results are presented in Fig.|6^, 
where we show P^o as a function of V{a) for the three different fits 
to the highly- variable hot pressure data presented by ISergis et"ai] 
( l2007l) . which they referred to (in order of increasing plasma f3) as 
the quiescent (blue), average (black) and disturbed (red) ring cur- 
rent. The shaded regions indicate the variability introduced by mod- 
ifying the subsolar standoff distance of th e magnetopause w hich 
modifies the parameters of the CAN model teunce et alj2007l) . By 
comparison between the calculations and the isotherms (solid) and 
adiabats (dashed) one can see that there is only a very narrow re- 
gion in 1/ ~ 12-16 Rs where the transport can be considered to 
be either isothermal or adiabatic. Inside 9 Rs Saturn's neutral OH 
cloud is a strong absorber of energetic particles and losses might 
reasonably account for the decrease in pressure at smaller values 
of L (and V{a)). In support of this, it is known th at the hot oxy- 
gen t emperature is approximately constant with L JPialvnas et al.l 
I2OO9I) and that the hot oxygen contributes the most to the hot pres- 
sure Pho- Hence, a reduction in pressure is related to a decrease in 
the number of hot oxygen ions per unit flux. At larger distances the 
pressure varies more steeply than PhoV{ay^^''^ suggesting a reduc- 
tion in pressure. T his may be related to the observed warping of the 
magnetic equator /Arridg e et alj2008 a^ which implies that particle 
pressures measured in the rotational equator will be smaller than 
at the magnetic equatoria l plane. Ene r getic particle pressures be- 
yond 20 Rs presented bv lSergis et 120091) support the fact that 
the pre ssure appears to be underestimated by the fits in lSergis et al.l 
( I2OO7I) . 

The product PhoV{a) was also determined as a function of p 
and shown in figure Fig.[6j3. As expected, PiioV{a) increases lin- 
early with L within, and just beyond Saturn's neutral cloud, due 
to the increasing flux tube volume and pressure and peaks near 
13-15 Rs be fore falling with increasing L. Using the pressures 
published by 1 Sergis et alj ( 12009) and the calculated flux tube vol- 
umes beyond around 16 Rs the value of PhoV{a) for L > 16 Rs 
is greater than ~ 5 • 10^ PamT~^. Nevertheless, the entirety of 
the PhoV{a) profiles reflects the strong variability in hot plasma 
f3 for Saturn's magnetosphere. It is important to note in this con- 
text that the differ ent ring current 'states' for plasma fitted by 
ISergis et alj j2007h represent the range of values of this parameter 
over several different orbits and magnetospheric configurations, as 
well as a wide range of radial distances and local times within each 
orbit. 

In view of the strong variability in this parameter and its gen- 
eral decline with decreasing distance inside ^ 10 Rs, we adopted 
a sim ple representation for its global behaviour, similar to that of 
ICau dal ( 1986). We composed a profile of Pm) by setting the prod- 
uct PhoV{a) to a constant value Kh beyond p — 8Rs; and by 
decreasing hot pressure linearly with decreasing p inside this dis- 
tance, according to the formula Pho{p) = Pho(8Rs) x (p/8). 
We then retrieved Pho values in our model's outer magnetosphere 
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Figure 6. Hot plasma pressure Pf^o as a function of unit flux tube vol- 
ume V{a) in SI units (panel a) and the product of hot plasma pressure 
and unit flux tube volume as a function of L (panel b) in Saturn's equa- 
torial plane. The coloured curv es employ the fitted pressure profiles de- 
termined by ISergis et alj j2007l) to the highly variable pressure data from 
MIMI. The colours represent the disturbed (red), average (black) and qui- 
escent (blue) ring current states, in the parlance of these authors. The solid 
coloured curves represent the profiles for a nominal magnetopause standoff 
distance of 25 Rs and the shaded regions represent the variability intro- 
duced into these profiles due to the changing upstream solar wind dynamic 
pressure affecting the global magnetic field configuration. In panel (a) the 
solid lines represent isotherms and the dashed lines ai'e adiabats. Particular 
L-shells are indicated by the solid dots and annotation. In panel (b) the hor- 
izontal dotted line shows the profile corresponding to the constant value of 
Kh = PV = 2 ■ 10*^ PamT-l (see text). This value may be compared 
with the data-derived curves, and we use this profile as a simple means of 
representing 'average conditions' within the highly-variable pressure distei- 
bution beyond p = 8 Rs . 
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beyond 8Rs through P^o ~ Kh/V{a). This form gives a more 
realistic response of the value of hot pressure to different configu- 
rations of the outer magnetosphere (from expanded to compressed) 
than would a single function of p alone. In addition, the param- 
eter Kh gives a compact representation of the 'level of activity' 
of the ring current in the model, and also reduces the number of 
free parameters. We intend to pursue a future parametric study of 
disc structure dependent on this parameter and magnetopause ra- 
dius. For the purpose of this introductory study, we set the value 
in our calculations (the scaling factor for 
normalised Kh is given in Table lBlt . This value according to Fig.|6] 
represents a ring current somewhat more 'disturbed' than the aver- 
age state. 

2.3.4 Plasma Rotation 

Eq. Isj for the plasma source function includes the scale length 
£, itself dependent on the angular velocity of the cold rotating 
plasma at each point along the equatorial plane. Under the steady- 
state assumption of the model, this angular velocity lu is constant 
along a magnetic field line, i.e. uj is expressible as a function of 
a alone. Thus it is observations of plasma ui which we seek in or- 
der to complete the equatorial boundary conditions for our model 
calculations. To construct a model for the azimuthal veloc ity va , 
of t he plasma we used data from st udies by Kane et al] ( 1200 8h 
and l'Wilson eTZI jlOOSl) . The study bv lKane et al.l ( |2008|) provided 
measurements of v,f, through analysis of ion velocity anisotropies, 
acquired in Saturn's outer magnetosphere by ion-neutral camera 
(INCA), a detector of MIMI. These data were acquired from the 
ion mode of the INCA instrument. As well as the estimated un- 
certain ty of ~ 20 per c ent in their individual velocity measure- 
ments, iKane et a\. Illooi's results also show considerable variabil- 
ity, around factors of two, within subsets of their measurements ac- 
quired near the same radial distance. This variability is attributable 
to the underlying set of spacecraft orbits sampling different local 
times and magnetospheric configurations (e.g. the INCA data used 
were obtained in the dawn sector for p < 25 Rs and in the mid- 
night sector outside this d istance). Following a different approach 
with a different data set, IWilson et alj (jlOOS*) determined by 
fitting drifting bi-Maxwellian velocity distributions to CAPS ion 
mass spectrometer data. They presented quadratic fits of (L) for 
the region between 5.5 and ~ 10 Rs. 

While a fully self-consistent model would include the influ- 
ence of magnetospheric configuration on the profile of and angu- 
lar velocity cu, we shall address this issue in a future study. For the 
present purpose, we use a profile of v^f, versus p obtai ned by fitting 
a sixth -order polynomial to points from the mo del of 'Wi lson et al.l 
( I2OO8I) and points taken from Figure 4 of I Kane et al. (2008}). Inside 
of 3.1414 Rs we assumed the plasma is in ideal corotation with 
an angular velocity of 1.638 • 10"* rad (a period of 10.65 h). 
For the purposes of fitting, we used a constant value of outside 
25 Rs equal to the average value of the outer magnetospheric ob- 
servations from Kane et al. (2()ol) and^cAndrews et al. (200^. 
We found that this approach produced a well-behaved fit in the 
outer magnetosphere without large 'oscillations', as well as good 
agreem ent in the inner magnetosphere with the data of lwilson et al.l 
( l2008n . For the calculations in this paper, we used the resulting 
polynomial fit to represent plasma angular velocity throughout the 
modelled magnetosphere. However we emphasise that fixing the 
value of to 169.25 kms~* beyond 25 Rs does not significantly 
alter the conclusions of our study. The and ui profile correspond- 
ing to this fit are illustrated in Fig. [7] For further comparison, we 
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Figure 7. Upper panel: A polynomial fit of order six for the azimuthal 
plasma velocity (thin solid curve) compai'ed with observations of the plasma 
azimuthal velocity in Saturn's magnetosphere. The squares are data from 
Kane et al. (2008) and were used for our fit. The thick, dark grey line is 
Wilson et al. (2008)'s empirical profile derived from their data, which were 
als o used for our fittin g. The light triangles were derived from Voyager data 
bv lRichardsor]< 19981) and are shown for compaiison only but were not used 
in our fit. Lower panel: The same comparison in the upper panel is shown 
with the azimuthal velocity transformed to angular velocity in the models 
and observations. 



include the Voyager velocity measurements by 'Richardson' ('199^ 
in the figure, but emphasise that we did not use these measurements 
to derive our fitted profiles. At distances smaller than ~ 10 Rs, the 
model curve agrees well with the Voyager data. Beyond this dis- 
tance, the model has values higher than the mean Voyager values, 
but is still consistent with the full range of these measurements. 

For the information of other modellers, we also present here 
the seven-element vector C of polynomial coefficients for the fitted 
plasma velocity profile. The following coefficients generate v^f, in 
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3 MAGNETODISC MODELS 

Having described our methods for incorporating equatorial obser- 
vations of plasma properties into the Caudalian model formalism, 
we now turn our attention to some example model outputs and how 
such calculations may be used to infer some important physical as- 
pects of magnetodisc structure at Saturn. We shall firstly consider 
some aspects of force balance in a disc formed under average solar 
wind pressure conditions and magnetopause size as observed in the 
Cassini era. 



3.1 Magnetodisc Structure for Average Magnetopause Size 

The probability distribution of magnetopau s e stan doff distance at 
Saturn was determined by lAchilleos et al.l jlOOSl) who surveyed 
magnetopause crossings of the Cassini spacecraft during 16 orbits 
between July 2004 and September 2005. The mean standoff dis- 
tance for this interval was found to be ~ 25 Rs . We thus adopt this 
value for our present work as an appropriate magnetopause radius 
for a nominal magnetodisc model representing average solar wind 
conditions at Saturn. The presence of the magnetopause boundary 
requires a corresponding contribution to the magnet i c pote ntial a 
from the currents flowing on that boundary. ICaudall repre- 
sented this magnetopause potential at Jupiter as the Euler function 
corresponding to a globally unifor m, southward- directed field Bs, 
referred to as the 'shielding field'. ICaudall ( Il986 ) chose the magni- 
tude of Bg by requiring that the magnetic flux due to the shield- 
ing field, integrated over the entire equatorial plane, be equal to 
a prescribed fraction ^ of the total magnetic flux exterior to the 
boundary due to the planetary plus disc sources. The addition of 
the shielding potential to the solution for a at each iteration thus 
'compresses' the flux tubes of the outer magnetosphere inwards 
from their 'boundary-free' configuration. 

For the magnetopause contribution in our axisymmetric mod- 
els, we adopted a similar approach to Caudal ( 1986); however, we 
determined our value of the uniform field Bs by performing day- 
side equatorial averages of th e empirical field mode ls described by 
lAlexeev & Belenkaval ( |2005|) : lAlexeev et alj ( |2006|) . which repre- 
sent contributions from both the magnetopause and magnetotail 
current sheets at Saturn. These two contributions are oppositely- 
directed (magnetopause field southward, magnetotail field north- 
ward). We computed our shielding field as a function of Rmp, 
using the following parameters to represent appro ximate condi- 
tions at Saturn, as required in the expressions of lAlexeev et al.l 
( I2OO6) : (i) planetary dipole orthogonal to the solar wind flow di- 
rection, (ii) radial distance Rt of the inner edge of the tail sheet 
equal to 0.7 Rmp, (iii) magnitude of field in the tail lobe given by 
Bl = $lRmp"^/ [|(1 + 2Rt/Rmp)], with open magnetic flux 
"l>i^ = 40 GWb. The resulting magnetopause contributions, before 
global averaging, showed variation by a factor ~ 2-3 between noon 
local time (strongest field) and the dawn / dusk meridian. The mag- 
netotail contribution showed similar relative variability but with the 
strongest fields situated at dawn / dusk. Thus in the full representa- 
tion there are local times where the two contributions add to zero. 
The uniform (dayside-averaged) shielding field used in our model 
is shown as a function of Rmp in Fig.[8l We show the contributions 
to the total shielding field from the magnetopause and tail currents. 
It is evident that for Saturn the magnetopause currents dominate 
the shielding field for the more compressed magnetosphere. For the 
more expanded configuration, the presence of tail currents signifi- 
cantly decreases the shielding field magnitude below its predicted 
values from magnetopause currents alone. 
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Figure 8. Dayside-averaged equatorial values of the shielding magnetic 
fields associated with Saturn's magnetotail current sheet ('Tail') and mag- 
netopause cun'ents ('MP'). The field value shown is the Z component, i.e. 
positive northwai'd. T he global fields used f or the averaging were calculated 
using the formulae of lAlexeev et al.l (2006) for the configuration where the 
planetary dipole is orthogonal to the upstream solar wind. 



Several output parameters associated with our average model 
(Rmp = 25 Rs) for the Kronian magnetodisc are depicted on a 
colour scale in the panels of Fig. |9] Fig. |9ja) shows contours of 
constant magnetic potential a, equivalent to magnetic field lines, 
for the vacuum dipole used to represent Saturn's internal field in 
our model. We may compare this geometry with the average mag- 
netodisc model in Fig.jgfb) which corresponds to magnetopause ra- 
dius Rmp = 25 Rs. The radial stretching of field lines compared to 
the dipole model becomes particularly pronounced beyond ~ 8 Rs . 
For example, the magnetic flux contained between the equatorial 
distances 6-10 Rs in the dipole field becomes spread out over a 
larger interval 8-18 Rs in the full magnetodisc solution. We shall 
compare the equatorial field profiles for these models later in this 
section. 

We now consider Fig. |9{c), which shows the distribution of 
total plasma pressure in the (p, Z) plane. The scale length I for the 
model ranges between 1-5 Rs through the magnetosphere, mono- 
tonically increasing with p. The pressure contours which attain sep- 
arations from the equatorial plane significantly larger than these 
scales are primarily due to the hot plasma pressure, which we have 
assumed to be uniformly distributed along field lines. One can also 
see the influence of the equatorial confinement of the cold popula- 
tion, by comparing individual contours with the field line shapes: 
the pressure contours tend to be more oblate. Fig. |9{d) shows the 
magnetic pressure distribution, along with contours of plasma /3, 
which clearly show the influence of the equatorial confinement of 
the cold population for /3 of the order unity or larger. The con- 
tours of magnetic pressure turn inwards towards the planet as they 
approach the equator. This is a consequence of force balance per- 
pendicular to the radial ly-stretched field lines ju s t outsi de the equa- 
torial plasma disc (e.g. iKivelson & Southwoo j |2005)). The main 
forces acting in this direction (which is approximately perpendicu- 
lar to the equator) are the plasma and magnetic pressure gradients. 
To maintain balance as the disc is approached, the corresponding 
increase in plasma pressure must be balanced by decrease in mag- 
netic pressure; hence the behaviour of the magnetic pressure con- 
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Figure 9. (a) Logarithm of magnetic potential a for vacuum dipole con- 
figuration, sliown on a colour scale as a function of cylindrical co-ordinates 
p and Z; (b) Logarithm of magnetic potential for an average magnetodisc 
field model for Saturn, with magnetopause radius 25 Rs and hot plasma in- 
dex Kf^ = 2 10^ Pa m T~ ^ (see text); (c) Distribution of plasma pressure 
within the model from (b); (d) Distribution of magnetic pressure within the 
model from (b), along with labelled contours of constant plasma /3 (thick 
white lines); (e) Total (plasma plus magnetic) pressure on a colour scale. 
The approximately vertical pressure levels near the equatorial plasma disc 
are a consequence of force balance perpendicular to the magnetic field (see 
text). 



tours. We thus expect total plasma plus magnetic pressure to be 
constant along the vertical direction near the disc. Fig.|9le) shows 
contours of this total pressure, and confirms that they follow direc- 
tions nearly perpendicular to the equator. 

We shall continue our present investigation of average plas- 
madisc structure at Saturn by considering the model's equatorial 
properties of magnetic field and force balance in Fig.[TO] The upper 
panel compares the equatorial profiles of magnetic field strength as- 
sociated with the planetary internal dipole, and with our full mag- 
netodisc solution for average magnetopause size. As for the sim- 
ple zeroth-order disc models ( i]2.2K the presence of the plasmadisc 
produces a total field profile somewhat weaker than the parent 
dipole for the regions closest to the planet, and stronger than dipole 
field beyond a characteristic transition distance. The middle panel 
of the figure shows equatorial profiles of magnetic pressure, cold 
plasma pressure and hot plasma pressure. We note that the mag- 
netic pressure exceeds that of the plasma for distances smaller than 
~ 10 Rs. The hot pressure is the dominant source for distances 
around ~ 15 Rs- The bottom panel of Fig. [Tol shows the equato- 



rial profile of the absolute value of the various volume forces. We 
emphasise here that we have used line thickness to indicate regions 
where radial forces are directed outward (thicker lines) or inward 
(thin lines). Over most of the model magnetosphere the curvature 
force is the principal, inward-directed (i.e. negative radial) force. 
The sum of all the radial forces in the equatorial plane has a magni- 
tude less than 0.2 percent of the local curvature force; this fraction 
thus provides some measure of the degree of accuracy with which 
the model can simulate perfect force balance. 

The bottom panel of Fig. [TO] also indicates which forces dom- 
inate the balance and determine disc structure in different regions 
of the equatorial magnetosphere. Throughout the magnetosphere, 
the magnetic curvature force is the strongest inward-directed force. 
For distances p > 15 Rs, centrifugal force is higher than plasma 
pressure gradients by factors up to five, and is therefore the sec- 
ond most important term in the disc's stress balance. Closer to the 
planet, for p ~ 6-12 Rs, centrifugal force and plasma pressure gra- 
dients are comparable in magnitude, and the disc's field structure is 
determined by both sources of radial stress in approximately equal 
measure. T hese calculations are in broad agreement with the con- 
clusions of 'Arridge e t al.l who used current sheet crossings 
to show that centrifugal and pressure gradient forces were approx- 
imately equal in magnitude at 20 Rs whereas the model shows the 
centrifugal forces slightly larger at about twice that of the pressure 
gradient forces. 

Our average Kronian disc model contains a hot plasma pres- 
sure distribution which is indicative of a 'mildly disturbed' ring 
current (see Fig. |6]l. We therefore would expect hot plasma pres- 
sure to play a more dominant role in magnetospheric force balance 
under conditions of so-c alled 'disturbed' ri ng current, as shown by 
the Cassini observations jSergis et al.l2007l) . We defer a detailed in- 
vestigation of this aspect to a future study, and concentrate here on 
modelling conditions characteristic of the mean level of observed 
hot pressure. 

Fig. [To] shows a small 'kink' in the magnetic force profiles 
around 8 Rs; this is due to the sharp linear decrease we have as- 
sumed for characterising the product of hot plasma pressure and 
unit flux tube volume ( ij2.3.3l l. The termination at this distance of 
the curve representing the outward-directed force due to hot plasma 
pressure confirms a sharp change in the sign of the hot pressure 
gradient; this feature in turn corresponds to the rapid decline with 
decreasing distance of the modelled hot plasma density. The kink 
feature is thus somewhat artificial, but does not affect the validity 
of the global features of our modelled force profiles. 

We now consider the inner magnetospheric region (/o < 6 Rs) 
depicted in Fig.[TO] Inside this distance, the cold plasma population 
density rapidly decreases (as also shown by the behaviour of the 
centrifugal force, which is proportional to cold plasma pressure). 
This magnetospheric region is then characterised by a relative ab- 
sence of plasma and a magnetic field dominated by the planetary 
dipole. 

We conclude this section on the average magnetodisc structure 
at Saturn by investigating the relationship between the previously- 
considered forces which act to create the magnetodisc geometry 
and the magnetospheric currents which flow in response to the pres- 
ence of those forces. On a microscopic scale, we expect the main 
azimuthal currents to arise from drift motions associated with: (i) 
finite plasma pressure (gradient and curvature drifts), (ii) centrifu- 
gal force associate d with plasma rotation (inertial current) (e.g. 
iBunce et al.l l l2007h ). The macroscopic formalism of the model al- 
lows an alternative identification of these currents from force bal- 
ance considerations, as follows. For our Saturn model, the relevant 
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Figure 10. Upper panel: Equatorial profiles of magnetic field strength for 
the planetary dipole alone and for the full magnetodisc solution (Rmp = 
25 Rg). Middle panel: equatorial model profiles of magnetic and plasma 
pressures (hot and cold). Bottom panel: normalised volume forces in the 
equatorial plane of the model, labelled according to line style. We show 
the absolute value of force. Thick lines indicate positive (outward) radial 
force, while thin lines show regions where the force is inwardly-directed 
(negative). 



scaling factor for is listed in Table IB ll In order to separate the 
contribution of a particular force to the current density, we sim- 
ply substitute its corresponding contribution to the source function 
Eq. ^ for the function g used in Eq. Q. 

Following this method, we calculated the various contribu- 
tions to azimuthal current in the equatorial plane of the average 
magnetodisc model. Profiles of positive (in the direction of plan- 
etary rotation) are shown on a logarithmic scale in the upper panel 
of Fig. [TT] It is clear that the force which is associated with the 
dominant contribution to the magnetospheric current depends on 
radial distance. For example, we note that there is a broad local 
maximum in the centrifugal inertial current centred at ~ 16 Rs. 
This feature corresponds to a similar local maximum in plasma an- 
gular velocity according to the model profile from Fig.l?] The lower 
panel of Fig. [TT] shows equatorial profiles of plasma j3 along with 
an equivalent 13 for the rotating disc plasma, computed as the ratio 
of the rotational kinetic energy density to the magnetic pressure. 
This rotational plasma /3 also peaks near ~ 16 Rg, thus indicating 
that the dominant term to the plasma source function (and therefore 
azimuthal current) in this region is the centrifugal force term. Rota- 
tional (5 then decreases for p > 16 due to the decline in cold plasma 
density. The corresponding effect on the current density profiles 
is a smaller ratio in the outer magnetosphere of the centrifugal to 
plasma current density. 

In the region p ~ 8-12 Rg, we see from Fig. II II that the 
current due to total plasma pressure gradient slightly exceeds the 
centrifugal current. The hot plasma current is an important factor 
here; the observed strong variability in hot plasma pressure at Sat- 
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Figure 11. Upper panel: Equatorial profiles of positive azimuthal current 
density taken from the Saturn magnetodisc model with Rmp = 25 Rg . The 
profiles are plotted on a logarithmic scale and are colour-coded according to 
the force with which they correspond in the plasmadisc's dynamical balance 
(hot / cold plasma pressure, centrifugal force). A function proportional to 
1/p is also shown in order to indicat e the form assumed for t he current den- 
sity in the CAN annular disc model IConnemev et alj fl98ll) : Lower panel: 
Equatorial profiles of plasma /3 taken from the Saturn magnetodisc model 
with Rmp = 25 Rg. Profiles are colour-coded according to the physical 
origin of the energy density used to compute the /3 ratio (hot / cold plasma 
pressure, rotational kinetic energy). The dashed and solid grey curves show 
the plasma /3 profi les fitted to observations of the hot plasma pressure at 
Saturn obtained bv lSergis et ai]j2007ll2009l) . 

urn ( ^2331ISergis et aljJioOThdKrimigis et aI.I (l2Q07')) implies that 
differing levels of ring current activity may plausibly increase the 
radial extent of this region where plasma pressure dominates mag- 
netospheric current, or even lead to its disappearance. Inside 8 Rg 
the hot plasma density sharply decreases and the corresponding de- 
crease in the associated current profile produces an inner region 
where centrifugal current is once more the major contribution. We 
shall defer a detailed investigation of the influence of hot plasma 
index Kh ( ^2.3.3l l on magnetospheric current profiles to a future 
study. For present purposes, we note that the calculations indicate 
it is expected to play a significant role in determining the extent 
of the region where hot plasma pressure is the major source of the 
azimuthal current density. 

Alongside the modelled plasma (5 in the lower panel of 
Fi g. I HI we also show obser ved values of hot plasma 13 presented 
by ISergis et a i] ( I200ll2009h (grey solid and dashed curves). The 
dashed grey curve s indicate fits to hot plasma observations by 
ISergis et al.1 l l2007h , which show the variation in hot plasma (3 be- 
tween average and disturbed ring current states (see also ii2.3.3t . 
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The solid grey curve was determined from a more recent fit to the 
median equatorial values of hot plasma (5 at Saturn determined by 
(computed over ~ 0.1 Rs intervals), who in- 
cluded the significant contribution (~ 50 percent) to hot pressure 
due to O"*" ions. If we compare this curve with the median profile 
from the earlier study (light grey dashed curve), we see that this 
inclusion has significantly increased the hot plasma /3 which would 
characterise an average state of the ring current. Given this result, 
and the large intrinsic variability in observed hot plasma pressures, 
we believe that our simplified model distribution of hot plasma j3 
(red curve) is in reasonable agreement with the expected global be- 
haviour of this parameter. For the region p > 20 Rs o ur modelled 
hot plasma 13 is in excess of the declining values of ISergis et al.l 
ll2007l.l2009h; however, the distant magnetospheric observations by 
iKrimigis et5] ( l2007h (see Fig.[6]( do show hot plasma /? which are 
consistent with our choice for Kh- Improved future determinations 
of plasma moments will no doubt enable us to further refine our 
plasma parametrisation , but for the present study we shall remain 
with the description given in i]2.3.3l 

Another important feature of the plasma /3 profiles in Fig. [TT] 
is the relatively uniform ratio of ~ 2-3 in the outer magneto- 
sphere (p > 15 Rs) between hot and cold plasma /3. Our cal- 
culations also show that in this region the length scale I for the 
cold disc plasma (Eq. ^) monotonically increases with distance 
between ~ 3-5 Rs . If we use these values in Eq. ^ for the tran- 
sition distance in a homogeneous plasma disc between pressure- 
and centrifugally-dominated regions, we obtain pT ~ 12-22 Rs. 
This range of transition distances is consistently smaller than the 
model magnetopause radius. The actual transition distance for the 
model appears to be situated at ~ 12 Rs, beyond which distance 
the centrifugal current persistently exceeds total plasma pressure 
current. The homogeneous disc predictions for pr thus suggest that 
the actual value of this distance is expected to lie within Saturn's 
magnetosphere, and the full magnetodisc model confirms that this 
is indeed the case. The formula in Eq. {8) therefore provides a rea- 
sonable means of estimating the order of magnitude of pr from ob- 
served and / or theoretical properties of the magnetospheric plasma. 

We now make a comparative investigation of general magne- 
todisc structure by comparing the profiles in Fig. [TT] for Saturn's 
plasmadisc with those shown in Fig. [12] for J upiter. The re sults 
in Fig. [T2] reproduce the model calculation by ICaudall l ll986h for 
a Jovian magnetodisc with magnetopause radius Rmp = 80 Rj . 
The most striking difference between the Jupiter model and Sat- 
urn model is the clear dominance of the Jovian outer magneto- 
sphere's equatorial current density by hot plasma pressure. The hot 
plasma current is the major contribution to total for distances 
beyond ~ 40 Rj. We also note a much stronger contrast between 
hot and cold plasma (3 for Jupiter compared to Saturn. While the ra- 
tio Ph/ Pc is an order of magnitude or more beyond ~ 40 Rj in the 
Jovian model, the same quantity is < 2 in the Kronian calculation. 
As a result, the current profiles due to hot and cold plasma pres- 
sure gradients show generally comparable values at Saturn, while 
at Jupiter the cold plasma current is an order of magnitude or more 
weaker compared to that of the hot plasma. 

These results indicate that the much more expanded magne- 
tosphere of Jupiter develops an outer region beyond ~ 40 Rj, 
where the cold plasma's angular velocity and density decline at 
a rate sufficiently rapid to produce a plasma whose main energy 
content arises from the thermal motions of the hot particle popu- 
lation. Near ~ 27 Rj in the Jovian model, the rotational plasma 
P exceeds the hot plasma /3 and the centrifugal current becomes 
comparable with the hot plasma current. This is qualitatively sim- 
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Figure 12. Upper panel: Equatorial profiles of positive azimuthal current 
density taken from the Jupiter magnetodisc model with Rmp = 80 Rj 
(reproduction of the calculation by Caudal 11986)). The profiles are plotted 
on a logai'ithmic scale, and ai'e colour-coded according to the force with 
which they correspond in the plasmadisc's dynamical balance (hot / cold 
plasma pressure, centrifugal force). A function proportional to 1/p is also 
shown in order to indicate the form assumed for t he current density in the 
CAN annular disc model IConnemev et alj il98ll) : Lower panel: Equato- 
rial profiles of plasma (3 taken from the Jupiter magnetodisc model with 
Rmp = 80 Rj. Profiles are colour-coded according to the physical ori- 
gin of the energy density used to compute the (3 ratio (hot / cold plasma 
pressure, rotational kinetic energy). 



ilar to the corresponding behaviour near ~ 16 Rs in the Kronian 
model. If we repeat the exercise of computing the transition dis- 
tance for the values of plasma P and scale length from the Jovian 
model, we obtain values of pr in excess of Rmp (the values for 
length scale are I = 5-40 Rj, increasing with distance). This in- 
dicates that the centrifugal current at Jupiter should never exceed 
the hot plasma current in the outer magnetosphere, according to the 
simple homogeneous disc model. The full magnetodisc model we 
have presented for Jupiter confirms this prediction, showing a hot 
plasma-dominated magnetospheric current beyond 40 Rj. 

We now consider the relative magnitudes of the magneto- 
spheric current at Jupiter and Saturn predicted by the models. Both 
Fig. II Hand Fig. 1121 show normalised current densities, expressed 
using scale factors of 280nAm~^ (Saturn) and 4800 nAm~^ 
(Jupiter) (see Table ISTt . Although the absolute value of the scale 
current at Jupiter is much higher because of that planet's stronger 
internal field, we note something interesting when we compare the 
normalised current densities at both planets within the same dis- 
tance range of < 25 planetary radii: the values of normalised at 
Saturn over 5-16 planetary radii exceed those at Jupiter by factors 
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of ^ 5. Since the distance scale is similar for botli models, we con- 
clude that this feature is an indication that Saturn's ring current pro- 
duces a stronger relative perturbation to the planet's internal dipole 
within this distance range. Interestingly, Vasyliunas ( 2008) arrived 
at a similar conclusion by considering the plasma outflows near the 
orbital distances of lo and Enceladus (~ 6 and ~ 4 planetary radii, 
respectively) and demonstrating that these flows would be expected 
to produce a stronger relative distortion of the planetary dipole for 
Saturn. 

Our model calculations also show a spatial profile of total 
in the outer magnetosphere, for both Jupiter and Saturn, which falls 
off more steeply with radial distance p than the 1/p dependence 
used by the CAN current disc model. This is an important point 
of comparison, as it indicates that an outer plasmadisc structure 
obeying radial stress balance has a characteristic spatial gradient in 
current density which is significantly different to that usually as- 
sumed in ring current modelling studies. Despite this difference, 
however, both the Caudalian and CAN disc models are suitable for 
reproducing the larger-scale observed structures in the magnetodisc 
field, as we shall see in the following sections. The main advantage 
of the Caudalian disc is that it also provides realistic spatial pro- 
files of current and radial force arising from self-consistent global 
distributions of plasma. 



Saturn Magnetodisc 
Compressed (Rn^p=18R ) Expanded (R|yp=30R ) 




3.2 Response of Magnetodisc to Solar Wind Pressure 

In this section, we parametrise the effect of solar wind dynamic 
pressure by varying the magnetopause radius Rmp in our model 
calculations. In Fig. [T3] we present model outputs calculated for 
two configurations. The first corresponds to strongly compressed 
conditions for the Kronian magnetosphere with Rmp = 18Rs, 
and the second is for a value Rmp = 30 Rs which is typical of the 
most expanded magnetospheric structures observed in the Cassini 
era ( Achilleos et al. 2008; Arridge et al. 2006). We emphasise that 
the plasma parameters of temperature, angular velocity, flux tube 
content and hot plasma index are identical in the two models. The 
final solution for the magnetic field within each model will change 
the mapping between these last two parameters and local quantities, 
such as number density and pressure, according to the frozen-in 
condition. 

We commence our comparison of the compressed and ex- 
panded magnetodisc structures by considering the top panels of 
Fig. [T3] which show contours of constant magnetic potential a, 
equivalent to field line shapes. The region of strongly-radial field 
near the equatorial plane, as seen in the average model ( ij3.1| >, is 
also present in the expanded disc, particularly in the range p ~ 
15-20 Rs. The compressed magnetodisc, on the other hand, dis- 
plays field line shapes which are far less radially 'stretched' and 
which more closely resemble the g eometry of a pure dipole (see 
Fig.|9}. A similar result was found bv lBunce et al.l ( l2008h who mod- 
elled the ring current for various magnetospheric configurations as 
revealed by Cassini magnetometer data from a selection of orbits. 
The colour scale of the upper panels in Fig. [T3] indicates that both 
compressed and expanded models have similar levels of magnetic 
flux threading their entire equatorial planes; we therefore expect 
higher field strengths to be present in the compressed disc. The 
middle panels confirm that this is the case. Equatorial profiles of 
total magnetic field strength relative to that of the planetary dipole 
are shown as a function of p. Beyond ~ 5Rs, the compressed 
disc model has a persistently stronger magnetic field than the ex- 
panded one. Around ~ 15 Rs, for example, the compressed field 



Figure 13. The left and right columns of plots con'espond respectively 
to Saturn disc models calculated for compressed (Rmp = 18 Rs) and 
expanded (Rmp = 25 Rs) configurations. Top panels: The logarithm of 
magnetic potential a is plotted on a colour scale for the labelled configura- 
tions. Middle panels: The equatorial ratio of total to dipole magnetic field 
strength is plotted for both magnetodisc configurations. The increased field 
strength of the compressed magnetodisc is appai'ent. Bottom panels: Equa- 
torial profiles of the absolute value of normalised volume forces for the 
compressed and expanded models, labelled according to line colour Line 
style is used to indicate the direction of the radial forces, with solid lines 
indicating outward force and dashed lines indicating inward force. 



has reached a magnitude twice as large as the expanded configura- 
tion. 

This behaviour of the field strength and geometry under 
strongly-compressed conditions has important consequences for 
the ensuing magnetic forces which operate within the plasmadisc. 
In the bottom panels of Fig. [13] we plot equatorial profiles of the 
volume forces due to plasma pressure gradients, magnetic pres- 
sure gradient, magnetic curvature and centrifugal force. The plots 
show that magnetic curvature is the principal, radially-inward force 
for both disc configurations. Closer inspections of the two curva- 
ture force profiles reveals a remarkable feature; the compressed 
model shows a stronger curvature force beyond ~ 8Rs, whose 
ratio with respect to the expanded disc attains a maximum of ~ 1.4 
at /9 ~ 15-17 Rs. The compressed model is able to maintain a 
stronger curvature force via higher magnetic field strength, despite 
the increased radius of curvature of the local field line. We also 
show plots of the total magnetic force J x B for both models (sum 
of curvature force and magnetic pressure gradient). A comparison 
of the two sets of curves reveals that, beyond ~ 8 Rs, the magnetic 
pressure gradient in the compressed disc is larger relative to the 
curvature force than in the expanded case. This beh a viour i s qual- 
itatively consistent with the study by lArridge et all ( l2008bh men- 
tioned in ^ which showed that the dayside magnetospheric field 
at Saturn only becomes significantly 'disc-like' under conditions of 
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low solar wind dynamic pressure (Rmp > 23 Rs). This aspe ct is 
also in accordance with the conclusions of lBunce et alj ( |2008|) . 

Within the range of radial distances 1 < p < 18Rs covered 
by the compressed model's equatorial plane, there are also signifi- 
cant differences in magnetic pressure gradient and centrifugal force 
with respect to the expanded model. Firstly, the magnetic pressure 
within this distance range falls off with distance more gradually in 
the compressed disc. For both configurations, power-law fits to the 
magnetic pressure, Pmag oc p~^^, were obtained for the interval 
10 < p < 15 Rs. The resulting indices were x = 2.80 ± 0.14 
(compressed) and x = 3.27 ± 0.10 (expanded), revealing that the 
expanded model field falls off slightly more rapidly than a pure 
dipole X = 3 in this region. However, a similar fit to the apparently 
more uniform part of the expanded field strength profile in the more 
distant magnetosphere 20 < p < 25 Rs yielded x = 1-12 ± 0.08. 
These results indicate that the compressed Kronian outer magneto- 
sphere is likely to be characterised by field strength gradient sim- 
ilar to that of a dipole, while a more expanded configuration may 
be expected to exhibit a field with a more gradual decline, asso- 
ciated with values of the index x in the range 1-3. This predicted 
behaviour of the magnetospheric field suggests that observational 
studies of the relationship between magnetopause standoff distance 
and solar wind pressure may benefit from the assumption of a 
field strength index x which varies with Rmp , rathe r than the usu- 
ally as s umed fixed value (e .g. lAchilleos et alj ( l2008l) : lArridge et al.l 
( l2006l) : ISlavinetalJjl985l) ). 

If we now turn our attention to the centrifugal force profiles 
in Fig.[T3] a detailed inspection reveals that the compressed model 
exhibits a centrifugal force consistently stronger than that of the 
expanded disc for ~ 8 < p < 18Rs, with the ratio of the two 
increasing monotonically to a value of ~ 2. This is a consequence 
of the higher cold plasma densities in the compressed model (at 
a given p, the ratio of centrifugal force between the two config- 
urations is equivalent to the ratio of cold plasma density). In the 
region ~8<p< 15 Rs, the plasma pressure gradients in the 
two models differ by less than 10 per cent. The increased cen- 
trifugal force of the compressed disc is thus balanced by an in- 
creased magnetic force (difference between magnetic curvature in- 
ward and magnetic pressure gradient outward). Interestingly, the 
region ~ 15 < p < 18 Rs near the compressed magnetopause 
is characterised by a change in sign of the magnetic pressure gra- 
dient, which is required to maintain balance due to the curvature 
force decreasing more rapidly than the centrifugal force. 

The region ^ 15 < p < 20 Rs for the expanded magnetodisc 
has force balance mainly determined by magnetic curvature and 
centrifugal effects as shown by Fig.[T3] The same region, however, 
has a broad local maximum in centrifugal force nearly coincident 
with a minimum in magnetic pressure gradient. These features arise 
because of the corresponding local maximum of plasma angular ve- 
locity in the same region (Fig.|7}, and the field geometry producing 
a relatively uniform region of field strength. In the more distant 
magnetosphere near ^ 23-27 Rs the hot plasma pressure gradient 
becomes equal in importance to centrifugal force in maintaining 
force balance due to the declining density and angular velocity of 
the cold disc plasma. As for the case of the compressed disc, the 
magnetic pressure near the magnetopause in the expanded model 
begins to increase with distance in order to maintain force balance 
near the boundary. 



3.3 Comparison of Model to Magnetic Field Observations 

A comprehensive comparison of our magnetodisc model for Saturn 
with the vast field and plasma datasets from the Cassini spacecraft 
is beyond the scope of this paper. For the sake of a preliminary 
assessment of how well the model may be applied to spacecraft ob- 
servations, we shall present in this section a comparison between 
the magnetodisc model field and magnetometer observations from 
the Cassini spacecraft from two quite different orbits. The first is 
the Revolution 3 (Rev 3) orbit lying entirely within Saturn's rota- 
tional equator, from the early part of the mission (February 2005) 
and the second is the highly-inclined Rev 40 orbit from March 2007 
which sampled the entire vertical structure of the disc. 

Fig. [14] depicts the information relevant for our comparison 
based on Cassini Rev 3, covering a period of approximately four 
days in February 2005. The time axis is labelled in days since 
the beginning of Day of Year 44, or February 13. The bottom 
panel shows the position of Cassini as a function of time using the 
colour-coded p and Z cylindrical coordinates as well as the Sat- 
urn local time (SLT) in decimal hours. We see that this inbound 
pass of the orbit sampled the magnetosphere at radial distances 
p ~ 4-31 Rs, the largest distance in this range corresponding 
to the indicated magnetopause crossing. The orbital segment for 
which p > 8Rs was situated at near-noon local times between 
■~ 10 and 14 h (SLT). The entire orbit was also situated within or 
very close to the equatorial plane Z — 0. This region of space is 
thus appropriate for analysis with our model, which is represen- 
tative of dayside conditions at Saturn and which uses a simplified 
formulation for the field due to magnetopause currents, based on an 
empirical model of the dayside equatorial field due to this source 

(313. 

The upper panels of Fig. [14] show the Bp and Bz compo- 
nents of the magnetic field (black curve) observed by Cassini dur- 
ing the relevant time interval, in nT units, from which we have 
subtracted the components of the internal field model described 
by Dougherty et al. (2005). The plotted data thus represents the 
magnetic field due to the external sources of the current disc and 
magnetospheric boundaries. To analyse the Bz observations, we 
chose two models. The first is the non-homogeneous part of our 
magnetodisc model (i.e. the total field model minus the planetary 
dipole term) with hot plasma index Kh ~ 2- 10*' Pam (repre- 
senting approximately average ring current activity at Saturn) and 
appropriate magnetopause radius Rmp = 30 Rs. The value for 
Rmp based on th e magnetopau s e cros sing location and the magne- 
topause model of lArridge et al.l l|2005) is 28 Rs — the use of either 
value did not significantly change the results. The Bz values shown 
by the blue curve were obtained by linear interpolation of the field 
values from our 2D model grid onto the spacecraft trajectory. We 
also sho w, using grey curve s, the predictions from the CAN model 
used by iBunce et alj (|2007) to analyse these data. The parameters 
for this model are the azimuthal scale current per unit radial length 
7o, the inner and outer edges of the annular model disc (a and b] 
and the disc half width in the Z direction (D). The parameter values 
we chose were the following, as determined by Bunce et al. (2007]): 
^loIo = 53.3 nT, a = 7 Rs, & = 20 Rs, -D = 2.5 Rs. 

The thin blue and grey curves show the corresponding contri- 
butions from the magnetopause shielding field to the different disc 
models. On the scale of the plot, the uniform shielding field of our 
model (see Fig.[8]l has a barely discernible magnitude of 0.09 nT. 
As for our Caudalian model, the shielding field for the CAN model 
was assumed to lie entirely in the Z direction, but was computed 
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Figure 14. The top two panels show a comparison between modified 
Cassini magnetometer data from the Rev 3 orbit ('MAG') and models 
which employ the Connemey ('CAN') and Caudalian ('MDISC') discs, 
with parameters as described in the text. The thin coloured curves in the top 
panel represent the small contributions in the models from magnetopause 
and tail cun'ents. Vertical (Bz) and radial (Bp) field components are shown 
as a function of time. The da ta are hourly av erages and have had the inter- 
nal field model described bv Doughertv et ai] f2005) subtracted. The middle 
panel uses a magnetic equator for the models which is displaced by 2 Rs 
north of the planet's rotational equator (without such a displacement, the 
models' predicted values for Bp would be identically zero for this equa- 
torial orbit). The bottom panel shows the spacecraft position as a function 
of time (cylindrical radial (p) and vertical (Z) distance, SLT). The vertical 
dashed lines indicate, from left to right, the positions of the last inbound 
magnetopause crossing, the outer edge of the CAN model disc (20 Rg) and 
the inner ed ge of the sa r ne mo del (7Rs). The CAN disc parameters were 
taken from lBunce et alj i2007l) , who fit a non-displaced model to observa- 
tions of Bz ■ 

using the following formula from lBunce et alj jlOOTh : 

= xT^, ■ ^2^^ 

This expression describes a shielding field which changes lin- 
early with X, the spatial coordinate associated with the axis which 
lies along the intersection of the equatorial plane and the noon- 
midnight meridian {X positive towards the Sun). The parame- 
ter values Bi = OnT and B2 ~ — l.llnT were chosen to 
fit the observed Bz values at the position of the magnetopause 
(Xi = 23.36 Rs) and the nightstde location with the minimum 
value of X (X2 = — 6.17Rs). It is evident that the shielding field 



for both models makes only a very minor contribution to the total 
predicted field except, for the CAN model, in the region adjacent 
to the magnetopause. 

Structure at a variety of timescales is evident in the observa- 
tions. The global nature of the disc models implies that they are 
suitable for analysing the largest scales, of the order half a day in 
time or a few planetary radii in p. If we firstly consider the Bz field, 
both models reasonably reproduce the overall trend seen in the ob- 
servations. Near the location of the outer edge of the CAN model 
(vertical line at ~ 2.3 days), we see that this model predicts a rela- 
tively sharp minimum in Bz due to the truncated nature of its cur- 
rent disc. The Caudalian disc with its extended current sheet makes 
a smoother transition in Bz through this region in better agree- 
ment with the observations. On the other hand, near the location 
of the inner edge of the CAN model (vertical line at ~ 3.7 days) 
the local peak in Bz displayed by this model fits the data more 
closely than the Caudalian disc, which rises to values more than 
twice that of the data within this inner region. This suggests a need 
for more accurate plasma inputs in the region p < 5Rs of our 
model, as discussed in i]2.3.1| and ij2.3.2| The local peak in the Bz 
data near ~ 4 days is most likely a signature of the 'camshaft' field 
at Saturn. This is a quasi-periodic modulation seen in the magnetic 
field whose physical origin remains to be unambiguously identified, 
but appears to be linked with field-aligned, azimuthally-modulated 
magnetosph eric currents flowing between the ionosphere and pla s- 
madisc (e.g. lSouthwood & KivelsoiJ ( l2007h : |Provan et all (12009')). 
We shall return to this aspect when we consider the high-latitude 
observations. 

Both disc models fail to agree with the Bp data in Fig. [14] 
This is because they have a north-south hemispherical symmetry, 
which by definition requires Bp — within the equatorial plane. 
The fact that Cassini observes a significantly non-zero Bp in Rev 3 
and many other equatorial orbits has been suggested to be the result 
of a non-planar plasmadisc structure; in particular, the bowl-shaped 
current sheet model explored by Arridge et al. (2008a) provides an 
explanation for these observations. Such a sheet morphology would 
be expected to arise in a magnetosphere where the planetary dipole 
is significantly non-orthogonal with respect to the upstream solar 
wind flow direction, as was the case during Rev 3, where the an- 
gle between these two directions was ~ 70 ° (the northern mag- 
netic pole being tilted away from the Sun). As a first, albeit crude, 
approximation to the ensuing field geometry, the model Bp val- 
ues shown are those corresponding to a model plasmadisc which 
has been displaced by a distance of 3 Rs north of Saturn's ro- 
tational equator. This displaceme nt is consistent with t he predic- 
tions of the current sheet model bv lArridge et"ai] ( l2008ah for a dis- 
tance p = 25 Rs and a subsolar latitude — 23 °, appropriate for the 
time of the Rev 3 orbit. The Bp field of our model in the region 
p — 20-30 Rs changed by < 0.5 nT when we changed the dis- 
placement by 1 Rs. We shall address the variability of the current 
sheet displacement with p in a future study. For present purposes, 
we use 2 Rs as a representative displacement for this outer magne- 
tospheric region. 

In the region between the magnetopause boundary and the 
neighbourhood of the outer edge of the CAN disc, both models 
and data are in reasonable agreement, confirming that a displaced 
planar disc is a useful representation of the local effects of the more 
realistic bowl-like shape. As we proceed closer to the planet along 
the spacecraft orbit towards the inner CAN disc edge, the observed 
Bp decreases in magnitude, consistent with the magnetic equator 
of the plasmadisc becoming aligned with the rotational equator; 
the displaced model values, unsurprisingly, do not fit the data in 
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this region. We see a local peak in the observed Bp near 4 Rs cor- 
responding to the similar feature in Bz- We noted that this peak 
in Bz changed by ~ 5 nT if we u sed a different internal field 
model for subtraction teurton et alT20 09). Thus we cannot provide 
a definitive explanation for this feature without further examination 
of the internal field models used in the near-planet region. 

In the panels of Fig. [Ts] we plot magnetic field components 
and spacecraft position as a function of time using the same scheme 
and conventions as Fig. [14] The time interval in question covers 
about nine days from the beginning of March 21, 2007 which cor- 
respond to the closest approach to Saturn and outbound segment 
of the Rev 40 orbit. The Z co-ordinate trace in spacecraft position 
shows that during this time Cassini probed regions up to 15 Rs 
from the equatorial plane, and the spacecraft latitude reached mag- 
nitudes of ~ 60°. In addition, at the time intervals near 6.5 and 
10 days, the spacecraft traversed the full extent in Z through the 
current sheet, a structure with typical vertical length scales of a few 
Rs. From the discussion in i]3.1| typical length scales along the p 
co-ordinate for the cold disc plasma are ~ 1-5 Rs , which there- 
fore provide an upper bound for the length scale along Z. Rev 40 
thus provides a very different view of the magnetosphere compared 
to the equatorial pass of Rev 3 and hence a good means of further 
testing the suitability of the disc models for magnetic analyses. 

We shall consider firstly the Bz data and model predictions 
in the top panel of Fig. [15] As for Rev 3, structure in the magnetic 
field on a variety of time scales is seen; in particular, the quasi- 
periodic (~ 10.75 hr) camshaft signal in Bz is clearly evident with 
typical amplitudes of the order 1 nT. We shall return to this fea- 
ture presently after discussing the larger-scale features in the field 
profile. The Caudalian and CAN disc models used for comparison 
purposes are shown as thick blue and grey curves. We chose the fol- 
lowing parameters for the CAN disc, obtained from a least-squares 
fit to the combination of both observed field components as dis- 
played in the figure: polo = 40 nT, a = 6.6 Rs, 6 = 18.6 Rs, 
D = 3.2 Rs. The thin grey curve shows the shielding field pro- 
file used in the CAN model, computed using the following param- 
eters for Eq. ((21): -Bi = 0.5 nT, B2 = 1.5 nT, Xi = 15 Rs, 
X2 ~ — 14Rs. The thin blue curve shows a modified version of 
the Caudalian disc, which we describe in more detail later in this 
section. We note that the n ightside subset of these data was mod- 
elled bv lKellett et all ( l2009h using the CAN formulation. These au- 
thors used a substantially thinner current disc (D = 0.4 Rs) and a 
correspondingly more intense current parameter (polo — 338 nT) 
to optimally fit the nightside field. 

The location of the outbound magnetopause crossing at 
(p, Z) — (26.8, 17.3)Rs corresponds to a subsolar standoff dis- 
ta nce Rmp = 24 Rs, u sing the axisymmetric magnetopause model 
of I Arridge et al.l j2006l) . However, we found that the Rronian disc 
model with a somewhat larger magnetopause radius Rmp = 30 Rs 
gave significantly better agreement with the observations, and it is 
the field profiles for this more expanded model which we have dis- 
played. Since the magnetopause crossing was at a relatively high 
altitude Z above the equator, this finding may indicate that the mag- 
netopause of Saturn exhibits polar flattening, although evidence for 
this requires further studies of similar mid- to high-latitude bound- 
ary crossings. 

For the interval spanning closest approach until the outbound 
magnetopause crossing, both the empirical CAN model and the 
physical Caudalian model reproduce the large-scale trend in Bz 
from the magnetometry. The Caudalian disc predicts a mean field 
about 1 nT weaker than that observed in the time interval after 
10 days. This feature suggests that a positive shielding field, sim- 
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Figure 15. Plots analogous to Fig. 1141 for the Cassini Rev 40 orbit. Pre- 
dictions for a rotating disc whose axis is tilted at 10° to that of the planet's 
rotation axis are also shown ('MDISC (tilt)'). The vertical dashed lines in- 
dicate the positions of the inner (6.6 Rs) and outer (18.6 Rs) radii of the 
CAN disc, whose parameters have been chosen to best fit the Bp data (see 
text). The time axis is truncated on the right at the first outbound magne- 
topause crossing. 

ilar to that employed for the CAN model, may be a more realis- 
tic choice for this pass than the uniform negative value —0.09 nT 
used in our model (see Fig. [8](. The camshaft signal in Bz is ob- 
served throughout this orbit. This field source, when added to the 
planetary dipole, has been suggested to be equivalent to that of 
a tilted, rotating dipole in th e outer magnetosphere p > 15 Rs 
( ISouthwood & Kivelsonl200 7^. In this picture, we would expect the 
magnetic equator of the outer Rronian plasmadisc to also be tilted 
relative to the rotational equator. As a preliminary exploration of 
this concept, we have plotted in Fig. [15] a thin blue curve show- 
ing the field profiles associated with a tilted, rotating plasmadisc. 
We computed these profiles by simply taking the original magne- 
todisc model and transforming it to a co-ordinate system where the 
model's axis of cylindrical symmetry is tilted at an angle of 10° 
with respect to the planet's rotation axis (the latter now being de- 
fined as the Z direction in accordance with the data). The orien- 
tation of the model symmetry axis was also allowed to vary with 
time such that the azimuthal angle of its projection onto the rota- 
tional equator corresponded to a regular rotation with a period of 
10.75 h. We note that the observed camshaft signal does not have 
a fixed period, but one which may drift in value by the order of a 
minute over time scales of the order of a year, as revealed by its 
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radio signature, e.g. lKurth et alj j2007ll2008h . We include the tilted 
disc calculations here simply to emphasise that such a model can- 
not be consistent with the observations in their entirety. Although 
the amplitude of tilted disc Bz fluctuations match the data on the 
outbound pass reasonably well for p > 10 Rs, they rapidly di- 
minish inside this region. By contrast, the observations show per- 
sistent field fluctuations throughout the orbit. These quasi-periodic 
fluctuations, i.e. the 'camshaft signal', the phase relations between 
the different compon ents, and their ori gin, have been the sub- 
ject o f much research I Espinos a & Doug hertv I2OOOI : ICowlev et al.l 
liood ISouthwood & KivelsonI 12007, : .Provan et al.l l2009h . These 
studies also highlight the difference between phase relations of the 
camshaft field components and those of a simple rotating, tilted 
disc. 



4 SUMMARY AND DISCUSSION 

We have introduced a new m odel for Saturn 's magnetodisc, based 
on an orig inal formalism by ICaudall ( Il986l) . The model formal- 
ism is based on the magnetostatic solution for an Euler potential 
consistent with global balance between plasma pressure gradient, 
centrifugal force and magnetic force (J x B) in a cylindrically 
symmetric system. Such an approach has the advantage of be- 
ing able to predict a self-consistent system of plasma properties, 
magnetospheric azimuthal currents and magnetic field. The equa- 
torial boundary condition for the model was provided by observa- 
tions from the Cassini spacecraft of hot and cold plasma pressure, 
and cold plasma density and temperature ( ^2.3\ and subsections). 
In order to formulate a model with realistic global behaviour, we 
adopted relatively simple functional forms for these physical pa- 
rameters. In this c ontext, t he un it flux tube volume concept was 
applied, following ICaudall ( Il986l) , in order to compute global dis- 
tributions of plasma whose density and pressure would respond ap- 
propriately to different magnetospheric radii, according to the be- 
haviour expected of a 'frozen-in' plasma. The empirical fits to hot 
plasma pressure, cold plasma composition and cold plasm a temper- 
ature bv lSergis et all J200ll2009h : lwilson et alj feOOSi) were em- 
ployed in order to achieve a reasonable representation for the model 
equator of average magnetospheri c conditions at Sa t urn. We also 
used a polynomial fit to the data bv lKane et alJ ( l2008h : IWilson et al] 
1 I2OO8 ) for plasma angular velocity. 

Before presenting the outputs from the full model, we consid- 
ered a simple toy model emphasising the largest angular scales of 
the magnetic potential for a homogeneous disc, characterised by 
constant plasma l3, constant plasma scale length £ and full corota- 
tion with the parent planet. This simple model was used to show 
the influence on magnetic field geometry to be expected when a 
rotating plasmadisc is added to a planetary dipole. In particular, 
hot plasma pressure generally inflates outer magnetospheric flux 
tubes to greater radial distances while the centrifugal confinement 
of the rotating cold plasma towards the equator gives rise to in- 
flated, relatively oblate field lines. This zeroth-order behaviour was 
consistent with our full model for Saturn's magnetodisc under aver- 
age internal (ring current activity) and external (solar wind) condi- 
tions. Our consideration of this baseline model, for which magne- 
topause radius Rmp ~ 25 Rs confirmed the radial 'stretching' of 
the unperturbed dipolar field lines as a result of the currents flow- 
ing mainly in the equatorial plasmadisc. The corresponding equa- 
torial field strength in the model shows a region where it falls below 
the unperturbed dipole value for distances ~ 5-15 Rs and exceeds 
the dipole value beyond this range, also displaying a comparatively 



more gradual decrease with radial distance in the outer magneto- 
sphere. All of these features are general characteristics which arise 
from adding the solenoid-like magnetic field of the disc (ring) cur- 
rent alone to the planetary dipole. 

Examination of the equatorial radial forces in the average Kro- 
nian disc model revealed that, for distances beyond ~ 15 Rs, the 
principal forces determining disc structure are the magnetic curva- 
ture and centrifugal forces. This characteristic distance is consis- 
tent with the simple formula for the 'transition distance' pr be- 
tween pressure- and centrifugal-dominated structure which arose 
from the zeroth-order disc treatment. This formula reveals the con- 
ditions under which pr is most likely to exceed the magnetopause 
radius, and consequently the plasmadisc can never have a force bal- 
ance dominated by centrifugal force. The relevant conditions are: 
(i) hot plasma /? is very high compared to the cold plasma (i.e. 
thermal energy is large compared to rotational kinetic energy), (ii) 
plasma angular velocity is adequately low, or (iii) for a given tem- 
perature of cold plasma, its density is small (such that the quantity 
I (ic becomes very large. 

Consideration of the equatorial, azimuthal current density 
in the average Kronian disc model revealed that centrifugal in- 
ertial current was the primary contribution for distances beyond 
^ 13 Rs- For the region 8-12 Rs, hot plasma and centrifugal 
current were predicted to be comparable. However, a further explo- 
ration of more active ring current states (future study) is likely to 
show that this interval in distance will expand, and the hot plasma 
current intensify, as the hot plasma index Kh is increased beyond 
values appropriate for average conditions at Saturn. 

A comparison of the azimuthal current profiles between the 
average Sat urn disc mode l and a reproduction of the Jovian mag- 
netodisc bv lCaudallll98^ (for which Rmp = 80 Rj ) was also re- 
vealing. In particular, the calculations confirmed that the increased 
Ph/Pc ratio at Jupiter endows this planet's magnetosphere with 
equatorial azimuthal current dominantly due to hot plasma pressure 
beyond a distance ~ 30 Rj. Within ~ 20-30 Rj, the centrifugal 
and hot plasma currents are of similar magnitude. 

The normalised (dimensionless) quantities adopted in our 
model enabled us to make an important comparison between the 
strength of the azimuthal current density at Saturn and Jupiter. The 
values of normalised equatorial at Saturn, according to our cal- 
culations, are expected to exceed those at Jupiter by factors of ~ 
5 within the distance range ~ 5-16 planetary radii. The implica- 
tion of this result is that, while the absolute strength of the Kronian 
currents is far weaker than their Jovian counterparts, the relative 
perturbation to Saturn's internal field in this distance range due to 
the disc current would exceed that at Jupiter. 

In i ]3.2l we examined the response of the Saturn disc model 
to conditions of compressed (Rmp = 18 Rs) and expanded 
(Rmp = 30 Rs) magnetospheric configuration. Both models gen- 
erally showed centrifugally-dominated force balance beyond ~ 
15 Rs, although the expanded disc shows comparable centrifugal 
force and hot plasma pressure gradient near ~ 25 Rs due to the de- 
cline in cold plasma angular velocity. Interestingly, the compressed 
disc is able to maintain a similar or stronger curvature force than 
the expanded one, despite having nearly dipole-shaped field lines: 
this property is a consequence of the higher field strengths attained 
in the compressed magnetospheric state. Around ~ 15 Rs, for ex- 
ample, the compressed model's equatorial field already reaches a 
magnitude twice as large as the expanded configuration. Consid- 
eration of the gradients in magnetic pressure in the compressed 
and expanded disc models indicated that the index x = ^b^' 
which characterises the relative change in field strength B per rel- 
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ative change in radial distance p, is likely to vary as a function of 
magnetopause standoff distance. This dependence recommends the 
corresponding use of a variable x in future observational studies of 
the response of Saturn's magnetopause boundary to changing solar 
wind conditions. 

Finally, in ij3.3l we presented model calculations for a Kro- 
nian disc model with Rmp = 30 Rs and average hot plasma in- 
dex. We compared our model predictions for vertical {Bz) and ra- 
dial {Bp) field components with magnetometer data from two of 
the orbits of the Cassini spacecraft's prime mission. We also pre- 
sen ted model calculations for appropriate CAN annular disc mod- 
els l lConnernev et al.iri98lh as part of this comparison. In general, 
both the Caudalian and CAN disc model were able to account for 
the general large-scale trends seen in the data-derived magnetic 
field due to the magnetodisc current alone. However, certain dis- 
crepancies between the models and the observations point to a need 
to use non-planar disc geometries in more detailed studies. The 
first of these discrepancies is the non-zero radial field observed by 
Cassini during the Rev 3 orbit considered herein, which cannot be 
explained by a disc field with north-sou t h hemi spheric symmetry. 
An observational study bv lArridge et al.l 120083") for many equato- 
rial orbits revealed that this is a repeatable signature, and is most 
likely associated with a bowl-shaped current sheet. 

The second important discrepancy between our model calcu- 
lations and the magnetometry is the presence of observed quasi- 
periodic fluctuations in the field, known as the camshaft signal. 
While our rotating, tilted disc model was able to qualitatively re- 
produce similar field fluctuations in the outer magnetosphere (p > 
15-20 Rs), it was clearly not capable of explaining the observed 
behaviour of the camshaft signal for regions closer to the planet. 
A general advantage of the Caudalian disc model in the context 
of data interpretation is that its more realistic plasma distribution 
yields smoother predicted changes in field orientation for space- 
craft passes through the current sheet. The CAN model predicts 
sharp peaks in field components during such transitions due to its 
assumption of an annular geometry with definitive boundaries for 
the current-carrying region; it also shows similar abrupt changes in 
field for the regions near the assumed iimer and outer edges where 
the current region is truncated. The discussion in i i3.1l revealed that 
the force balance used to derive the Caudalian disc structure also 
results in a fall-off in magnetospheric current density more rapid 
than the 1/p law assumed in the CAN model. 

The Caudalian magnetodisc for Saturn represents a useful first 
model for pursuing studies of the plasmadisc structure, azimuthal 
current and magnetospheric field, along the lines that we have pre- 
sented in this paper. While these initial studies have revealed some 
interesting features of disc structure and currents at Saturn, particu- 
larly when compared to the Jovian system, they also highlight some 
important future directions for work involving this model, such as 
the following. 

(i) Improved determinations of plasma moments should be in- 
corporated into the structure of the model, in order to provide more 
accurate depictions of the global plasma conditions. 

(ii) Investigation of the influence of hot plasma pressure on mag- 
netodisc structure. Our initial study has revealed that it plays a 
potentially important role in determining the general structure of 
the magnetodisc field and the extent of the magnetospheric region 
where the electric current density is dominantly determined by 
energetic particle motions, rather than the inertial current associ- 
ated with centrifugal force acting on the cold population. 

(iii) An extension of our preliminary study of solar wind influ- 



ence on disc structure to additional magnetopause radii and global 
characterisations of internal plasma energy and content. The self- 
consistent response of plasma angular velocity to magnetospheric 
compression could also play a potentially important role here. 

(iv) Further analyses of Cassini field and plasma data. The 
spacecraft has thus far completed more than 100 orbits of Saturn. 
Such a vast dataset will require much time to exploit. A suitable 
use for our model with regard to the field and particle data would 
be a modelling study of selected orbits using input plasma moments 
acquired during those orbits, rather than a 'global approximation' 
to these conditions. The model outputs would thus reflect condi- 
tions most appropriate for the orbits in question. Such calculations 
would be of use, for example, to teams who aim to derive magne- 
tospheric particle fluxes and current densities directly from in situ 
measurements. 
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APPENDIX A: SOLUTIONS FOR THE MAGNETODISC 
POTENTIAL 

This Appendix describes the derivation of the solution for the mag- 
netic pote ntial of an ax isy mmetric p l asma distribution given in ar- 
ticles bv lCaudall ( Il986h and lLackned i ll970h . The potential in ques- 
tion is denoted by a and is actually one of two Euler potentials 
from which the magnetic field B may be derived 

B = Vq X V/3, (Al) 

where a is a function of radial distance r and cosine of colatitude 
^ = cos 9. The function (3 depends only on azimuthal angle 0: 
(3 = acf) with a being the planetary radius. Note that many pairs of 
Euler potentials can be associated with a particular magnetic field, 
however this particular choice separates the azimuthal and merid- 
ional dependencies. In effect, these equations tell us that an indi- 
vidual magnetic field line can be thought of as the line of intersec- 
tion of a surface of constant a (which will resemble a 'doughnut- 
shaped' shell) and a plane of constant (3 (which is simply the merid- 
ional plane with azimuth 0). 

Let us now consider the di fferential equa tion for the merid- 
ional Euler potential a given bv lCaudall ( fT983) ( we shall use the 
dimensionless system of co-ordinates described in this paper): 

In Eq. l lA2b . the function g represents a source term describ- 
ing a distribution of external plasma and currents which must be 
specified a priori. Note that g requires knowledge of a as well — 
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the fun ction we are trying to solve for. lCaudail ( Il986h and lLackneJ 
solve this problem through an iterative approach. One starts 
with an initial 'guess' ao for the functional form of a. ao is then 
used to evaluate g, and then Eq. iA2\ is solved to give an updated 
solution ai. ai is then used in the next iteration to re-evaluate g 
and to update the solution again. The process is repeated until con- 
vergence: in practice, one usually stops when the maximum rela- 
tive difference between successive iterations falls below some user- 
defined tolerance. 

A reasonable first guess for a is the planetary dipole potential 



Qdip(r,^) = 



(A3) 



Odip is a homogeneous solution of Eq. l lA2b (i.e. a solution for 
the case where the source term is identically zero). However, it is 
not the only homogeneous solution. Homogeneous solutions are 
a good starting point for finding particular solutions (i.e. when 
the source term is non-zero). We may obtain the general form of 
the homogeneous solution by using the property of separability i.e. 
a{r, = oir{r)a0{fj) is the product of two single-variable func- 
tions as stated previously. Substituting this into Eq. ( IA2t we can 
show that 



d ct, 



ag dy? 



(A4) 



Now if we fix the value of r, we would expect the left-hand 
term in this equation (a function of r only) to be a constant. How- 
ever the equation tells us that this constant is independent of what- 
ever value of y we use to evaluate the right-hand term. It follows 
that the right-hand term, regardless of the value of p, must be a 
constant. A similar argument, keeping y fixed, reveals that the left- 
hand term, for all values of r, must also be equal to a constant. If the 
constant takes on special values, derivable from an integer n ^ 0, 
we may write 



d Qr 



ar dr^ 
1 — fi'^ d^ae 



71(71 + 3) +2, 
71(71 + 3) +2. 



(A5) 
(A6) 



It is easy to show that the radial part has a solution of the 



form 



(A7) 



where C; is a constant, and integer I must satisfy 1(1 — V) = 71(71 + 
3) + 2. Solving this quadratic, we see that I can take on the value 
n + 2 or — (n + 1). / = 71 + 2 corresponds to a positive power 
of r and a potential which monotonically increases with distance 
from the planet — this is not physical. Therefore, we choose I = 
— (n + 1) for the radial part of the function. 

To solve for the angular function ag we require knowledge of 
the Jacobi polynomials. The particular strand of these polynomials 
which are of use to us here are denoted P^'^ (jj.) (n is an integer ^ 
so we can associate a polynomial with each choice of n in equations 
|A5| and |A6t . The useful property of the polynomials (/i) is that 
the functions (1 — fi^) P^'^{fi) are actually solutions of Eq. iA6i . 
For a given choice of n our homogeneous solution would thus be 



aH{r,iJ.) = QrQe = C„r 



'("+1) 



(1-m')P„'''(m)- (A8) 



Since Eq. l lA4t is linear in (a,.ae), it follows that any linear 
combination of solutions of the above form is also a homogeneous 



solution. Without loss of generality, the complete homogeneous so- 
lution is thus 

oc 

a„{r, p) = (1 - 11^) J2 Cur-^"+'^ P^'\y) . (A9) 

n=0 

We have now found solutions for the homogeneous (source- 
free) version of Caudal's equation. But how do we use these to 
obtain a solution for the full differential Eq. l|3j which contains the 
source function gl We try a general solution obtained by multiply- 
ing each term in the series of the homogeneous solution by a purely 
radial function fn{r). This trial function thus takes the form 

a{r,y) = (1 - y') J2 Mr)r~^"+'^ P^-^(m) , (AlO) 

where we have absorbed the constant C'„ into the definition of 
fn{r). Now if we use this trial solution in the left-hand side of 
Eq. ([3), we obtain 

n=0 ^ ^ 



: g{r,n,OLi- 



(All) 



Here we have introduced the symbol i to emphasise that solv- 
ing this equation is part of an iterative process where the solution 
Qi is obtained from the previous one a^-i. Although the left-hand 
side of our equation retains the form of a series summation in the 
Jacobi polynomials, we can't progress much further without ad- 
dressing the right-hand side. This is where the Jacobi polynomials 
again prove useful. They are an orthogonal, complete set of func- 
tions, which means that any function of y can be expressed as a 
series expansion using Jacobi polynomials. Applying this to our 
function g for an arbitrary value of radial distance r, we can decom- 
pose the angular dependence of g into a sum over the polynomials 
as follows 

00 

5(r,/x,ai_i) = (l-M^)^p„(r)P„i'i(/i), (A12) 

with the expansion coefficients defined by the orthogonality condi- 
tion 



ffn(r) 



g{r,y)Pi'\y)dy, (A13) 



hr. = j {l~y'){P^''{p)rdy. (A14) 

We can now make use of the orthogonality of the polynomials 
to equate the 71-th terms of Eqs. l lAllljAT2] l. This gives 

dV„ 2(n + l) df„\ 



r-^"+^) = -5n(r), (A15) 



dr^ r dr 

or equivalently (multiplying both sides by 7-~("+^)) 



2(n + l) d.^ fn 

dr^ 



- 2(71 + l)r-(2"+«) ^^ = 



dr 



(A 16) 



We see that the left-hand side can be expressed as the deriva- 
tive of a product 



(A17) 



The left-hand side of this equation is readily integrable. But we see 
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that the general solution for the fn (r) functions will involve inte- 
grals of the source function g. What this means in practice is that 
we have to numerically integrate some kind of empirical or other 
functio n which is a fit to observed plasma distributions. ICaudall 
( Il986h 's work shows that the source function includes quantities 
such as plasma pressure, plasma temperature (assumed isotropic) 
and mean ion mass. We now finalise the integration towards a final 
solution. We start with Eq. l lA17t and rename the dummy variable 
for radial distance to u 
d 



du 



(n + 1) 



gniu). 



(A 18) 



We now integrate both sides over the range to r. is an 
inner boundary, similar to the planetary radius, which encloses the 
region where the field is purely a dipole field i.e. purely due to the 
planet's internal source. We adopt the boundary condition that /„ = 
at M = rc (i.e. the contributions to the potential from the plasma 
source disappear at the inner boundary) and = at u = Vc 
(there is a 'jump' in the potential gradient at the inner boundary 
u = rc supported by currents flowing on that surface). Performing 
this integration between u — Vc and u = r gives us 



where G(r) denotes the function 

G{r) = r u-^"+'^g„{u)du, 
dG 



-(n + 1) 



dr 



(A19) 



(A20) 



(A21) 



If we integrate Eq. JAlSt between the limits u — rc and u = 
oo we obtain the useful identity 



We can now integrate Eq. JA19t by parts using the boundary 
conditions G{rc) — and fnij'c) = to get 



1 



2(n + l) 



du- 



2n + 3 



G(r 



u"~^^ gn{u)du 



(A23) 



2n + 3 

If we now make use of Eq. l lA22t to eliminate the unknown 
fc, and perform the first integral, we obtain 

1 



Mr) 



2n+3 



2n + 3 

1 

2n + 3 



2n + 3 



G{r) 



G(oo)- 

u^^'^ gn{u)du 



(A24) 



Since, by definition G(oo) — G(r) = J°° gnU''^"^^^ du, we 
can combine the two terms with factor r^"'^^ and multiply both 
sides by r"'"^^' to get the following form for the full radial part 
of the solution 



/nWr 



1 



2n + 3 

2n+3 



n + 2 



u ^"^^^ g„{u)du+ 



n + 2 

u g„ 



u)du- 
u"''"'^^'' gn{u)du 



(A25) 



We have written the solution in this form so that it reflects 
the full radial part of the solution given in Eq. l lAlOll. Thi s radia l 
part of the full solution agrees with that given bv lCaudall JT98i) 



Table Bl. Scaling values between variables in physical units and their di- 
mensionless counterparts for both planets. 



Dimension 


Definition 


Units 


Saturn 


Jupiter 




Primary scales 






Length 


a 


km 


60280 


71492 


Magnetic Field 


Bo 


nT 


21160 


428000 




Derived scales 






Volume 


a? 


km'' 


21014 


4-101-* 


Magnetic Flux 


Boa^ 


GWb 


77 


2187 


Cun'ent density 


Bo/(a/io) 


nA m^^ 


280 


4800 


Pressure 




Pa 


0.00036 


0.146 


Energy density 


Bo /mo 


J m"'' 


0.00036 


0.146 


Hot plasma index Kf^ 




PamT-l 


109 


2 ■ IQio 



and lLackned ( Il97(]h . Their work shows that the integral multiplied 
by a factor rc""*"^ comes about by assuming a boundary condition 
for fn different from zero. However, Caudal points out that this 
extra integral corresponds to surface currents at r = rc and makes 
negligible contribution to the solution beyond a few planetary radii. 
In fact for his final calculations he omits it and relies on a more 
detailed internal field model. For the work described in this paper, 
we use a simple centred dipole representation of Saturn's field, with 
equatorial field strength as given in Table[T] 

The final solution consists of the homogeneous part (assumed 
to be the dipole or other appropriate potential) added to the partic- 
ular solution (non-zero source) whose radial and angular parts we 
have derived above. For completeness, we now give here the final 
solution for the magnetodisc potential 



a{r,fJ.) = h 



n+2 



2n + 3 

, -(n + 1) / n + 2 / N I 

+ r ^ ^ ' I u ^ g„{u)du 

roc 

„2n + 3 / + g^^^^du 



gn{u)u ^"'^^^ du 



(A26) 



This represents, in practice, a cumbersome calculation. The 
number of terms required in the polynomial series depends on how 
accurate a representation is needed for the source function (whether 
empirical or theoretical). Source functions characterised by larger 
angular scales require fewer polynomials in the expansion. For the 
work described in this paper, we used polynomial expansion up to 
degree n = 30. The corresponding latitudinal resolution captured 
by the polynomial of this degree is ~ 2.2° , corresponding to typical 
vertical resolutions at the equator in the range 0.2-1 Ks - To obtain 
final model outputs, we stopped iteration when the maximum rela- 
tive difference in the solution for the magnetic potential a between 
consecutive iterations became less than 0.5 percent. 



APPENDIX B: SCALING OF PHYSICAL QUANTITIES 

Table IbTI presents a summary of the scaling values for all dimen- 
sionless quantities. 



24 A^. Achilleos, P. Guio and C. S. Arridge 



APPENDIX C: 
LENGTH 



DERIVATION OF THE PLASMA SCALE 



This appendix describes the derivation of the plasma scale length 
I defined in Eq. l|6j. Let us first consider the force balance along 
a magnetic field line assuming that the ions and the electrons are 
subject to an ambipolar electric potential <1?||. In these conditions 
the force balance equation for the ions and the electrons are written 
respectively 



+ nrriiUj pcos ip — ne 

ds ds 

2 , d^H 

+ nnieUJ pcosif + ne- 



0, 



= 0. 



(CI) 
(C2) 



_^ _ ^ _ 

ds ' ' 

Here s represents the curvilinear coordinate along a magnetic field 
line (i.e. such that the potential q(s) is constant) and is oriented to- 
ward the equator (i.e. such that positive force means force acting to- 
ward the equator). The angle tp denotes the angle between the mag- 
netic field line (in direction of increasing s) and the radial direction 
(with unit vector ep). In these force balance equations we have im- 
plicitly assumed quasi-neutrality of the plasma (i.e. Ui = = n), 
and also assumed that ion and electron pressure along the magnetic 
field are equal (i.e. Py . = Py^ = Py = nfci,r||). This ass umption 
is con sistent with the definition for total pressure used bv lCaudall 
( Il986l) for a quasi-neutral plasma with the specific ion charge num- 
ber Z = 1. However, it is worth noting that this simplifying as- 
sumption would have to be relaxed for more realistic studies of the 
observed differences between ion and electron pressures. Subtract- 
ing Eq. ( |C2t from Eq. JCll ) we obtain the following expression for 
the ambipolar electric field By, which is directed along the mag- 
netic field 

1 2 

S|| = — ^~-o — to P cos (fi. (C3) 
ds 2 e 

We note that is negative which means that the ions tend to be 
'lifted' off the equator to a greater degree than they would be in the 
absence of ambipolar effects. Substituting back the expression for 
the ambipolar electric field E^^ into the ion force balance equation 
Eq. iCl\ , we obtain the following differential equation 

dP|| 1 P|| p cos (f 
~d7 ~ ~2kb^/{miLLi'^ 



(C4) 



Changing the curvilinear coordinate s into the cylindrical radial 
distance p (i.e. dp — ds cos ip along the field line) we obtain the 
separable differential equation 

^ = PjP (C5) 

And finally expressing the radial distance p in normalised units and 
integrating from p to the equatorial crossing po (with pressure P||q) 
we obtain the following analytic expression for ion and electron 
pressure along a magnetic field line as function of radial distance p 



2P 



(C6) 



where we recognise the pla sma scale lengt h £ as given by Eq. l|6ll. 
Thus this definition of £ in ICaudall il9m 's formalism implicitly 
represents the above simplified expression for the ambipolar elec- 
tric field in the plasma. For our model, the charge state Z = 1 
and the value of ^ is a factor of y/2 larger than that which would 
be derived for ions in the absence of the ambipolar electric field. 
For plasma ions with general charge number Z, this factor is 
^/ZTT. As quasi-neutrality is maintained, the electrons are also 



distributed with the same scale length £ (since they have negligi- 
ble mass compared to the ions). A more thorough treatment of the 
polarisation electric field in the magnetospheres of Jupiter and Sat- 
urn, along with its ef fects on plasma distributions, can be found in 
lMauriceetai] ( ll997h . 



